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PICOSECOND  LASER  PULSE  INTERACTIONS  WITH 
METALLIC  AND  SEMICONDUCTING  SURFACES 

Principal  Investigators 
Eric  Mazur  and  Frans  Spaepen 


Preface 

This  report  contains  the  results  of  a  two-year  study  of  the  interaction  of 
picosecond  laser  pulses  with  metallic  and  semiconducting  srufaces,  and 
consists  for  the  most  part  of  a  collection  of  published  papers. 

The  study'  has  concentrated  on  the  time -dependence  of  electronic  and 
structural  transitions  on  a  picosecond-to-nanosecond  timescale,  following 
irradiation  with  a  picosecond  pulse. 

The  main  findings  of  the  work  are: 

1.  Use  of  streak  camera  with  1.8  ps  resolution  allows  the  reflectivity  during 
the  melting  of  silicon  to  be  studied  with  improved  spatial  and  temporal 
resolution  following  a  single  pump  pulse. 

2.  Following  a  30  ps,  532  nm  pump  pulse  at  470  mJ/cm2,  streak  camera 
measurements  show  that  silicon  overheats  above  its  melting  temperature 
by  more  than  1000K. 

3.  Pump-probe  and  ellipsometry  measurements  on  graphite  show  a 
substantial  decrease  of  both  the  imaginary  and  real  parts  of  the  index  of 
refraction  when  the  pump  fluence  exceeds  the  melting  threshold, 
indicating  that  the  molten  graphite  is  less  metallic. 

4.  The  use  of  these  short  probing  and  excitation  times  in  the  graphite  studies 
avoids  artifacts  arising  from  evaporation. 

5.  Pump-probe  reflectivity  measurements  on  the  layered  compounds  !T-TiS2 
and  lT-TiSe2  show  a  transition  to  a  phase  with  lower  reflectivity,  indicating 
that  is  not  metallic. 


6.  An  analysis  of  the  stability  of  the  liquid  created  by  picosecond  laser  pulse 
irradiation  shows  that  on  this  time  scale: 

(i)  The  loss  of  mass  and  energy  by  surface  evaporation  is  negligible. 

(ii)  The  steady-state  nucleation  rate  of  vapor  bubbles  is  sufficiently  small 
to  allow  substantial  overheating  of  the  liquid  above  its  boiling  point. 

(iii)  Transient  effects  can  play  a  role  in  postponing  bubble  nucleation. 

(iv)  The  stability  limit  of  the  overheated  liquid  (spinodal  decomposition 
to  the  vapor)  can  be  approached  to  within  a  few  hundered  degrees, 
corresponding  to  overheats  of  several  thousand  degrees. 

(v)  The  large  amounts  of  mass  transport  observed  after  the  laser  pulse 
are  not  the  result  of  evaporation,  but  are  caused  by  hydrodynamic 
effects  resulting  from  the  recoil  pressure  induced  by  surface 
evaporation. 

7.  Pump-probe  measurements  of  the  time  elemental  metals  (Au,  Cu)  stay 
molten  give  crystallization  velocities  up  to  100  m/s.  These  can  only  be 
obtained  by  a  collision-limited  crystallization  mechanism. 

8.  Quenching  of  Nb-Si  following  laser  irradiation  provides  a  case  study  of 
the  kinetic  preference  for  crystal  formation.  It  shows,  for  a  sequence  of 
five  phases,  how  higher  isotropy  and  a  smaller  unit  cell  size  are  favored  at 
higher  quench  rates. 
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Abstract 


The  use  of  a  streak  camera  makes  it  possible  to  improve  the  study  of  the 
melting  of  silicon  during  picosecond  laser  annealing.  Its  time  resolution  of 
1.8  ps  enables  us  to  confirm  that  the  molten  silicon  is  heated  above  the 
melting  temperature.  It  also  provides  spatial  information  on  the  melting 

process. 


1.  Introduction 


Numerous  investigations  on  the  phase  transition  of  silicon  during  ultrafast  laser 
annealing  have  been  performed  in  recent  years.  It  has  been  established  that  the  silicon 
surface  melts  during  a  picosecond  laser  pulse.1  Because  liquid  silicon  is  a  metal,  the 
reflectivity  increases  on  melting.  This  has  indeed  been  observed  using  optical  pump-and- 
probe  techniques.1'3  Standard  picosecond  pump-and-probe  measurements,  however, 
have  some  inherent  drawbacks.  First,  they  cannot  resolve  reflectivity  changes  that  occur 
on  a  time  scale  of  a  few  picoseconds,  because  they  integrate  over  the  duration  of  the 
probe  pulse  (typically  20  ps  or  more).  Second,  they  determine  the  time  profile  of  the 
reflectivity  for  every  pump  fluence  in  a  step-wise  manner  by  varying  the  delay  between 
the  pump  and  the  probe  pulse.  This  introduces  a  large  amount  of  scatter  in  the  data 
points,  due  to  shot-to-shot  variations  in  the  pump  fluence,  and  requires  a  large  amount  of 
data  to  be  taken  for  every  time  profile  of  the  reflectivity.  Also,  they  provide  no  spatial 
information  on  the  melting  process. 

To  obtain  both  spatial  resolution  and  better  time  resolution,  and  to  measure  the 
time  profile  of  the  reflectivity  on  a  single-shot  basis,  we  use  a  streak  camera  with  a  time 
resolution  of  1.8  ps  to  detect  the  probe  pulse. 


2.  Experimental  setup 

The  experimental  setup  is  shown  in  Fig.  1.  The  frequency-doubled  output  (30 
ps,  532  nm)  of  a  mode-locked  Nd:YAG  laser  is  split  into  a  probe  and  a  pump  pulse.  The 
duration  of  the  probe  pulse  is  stretched  by  splitting  it  in  four,  delaying  the  four  resulting 
pulses  with  respect  to  each  other,  and  recombining  them  spatially  to  form  a  longer  probe 
pulse  of  120  ps  duration.  The  probe  pulse  then  images  an  area  around  the  100-pm  wide 
melting  area  onto  the  entrance  slit  of  a  streak  camera  (Hamamatsu  Photonics  Cl 587). 

The  entrance  slit  of  the  camera  is  split  into  two  parts.  The  larger  part  is  used  to 
image  the  probe  pulse,  the  smaller  part  to  image  a  fraction  of  the  pump  pulse.  The  latter 
serves  as  the  timing  reference. 

To  enhance  the  sensitivity  of  our  measurements,  the  probe  pulse  is  p-polarized 
and  the  probe  angle  of  incidence  is  65°,  close  to  Brewster’s  angle.  At  this  angle  the 
reflectivity  of  solid  silicon  is  small  (about  10%),  leading  to  an  increase  in  reflectivity  on 
melting  by  a  factor  of  8. 


3.  Results  and  Discussion 

The  streak  images  of  two  measurements  are  shown  in  Figs.  2a  and  2b.  The  time 
axis  is  displayed  vertically,  with  time  increasing  from  top  to  bottom.  The  full  height  of 
the  image  corresponds  to  190  ps.  The  horizontal  axis  reflects  the  spatial  profile.  The 
pump  pulse  is  shown  on  the  left  side  and  has  a  duration  of  about  30  ps.  The  stretched 
probe  pulse  covers  most  of  the  image.  The  bright  part  is  where  the  silicon  surface  is 
melting.  Spatially  it  reflects  the  Gaussian  intensity  profile  of  the  pump  pulse.  The  time 
profile  allows  one  to  study  the  melting  dynamics. 

Both  measurements  were  performed  at  the  same  laser  fluence.  In  the  case  shown 
in  Fig.  2b,  some  surface  irregularity  leads  to  the  formation  of  a  surface  plasma  in  the 
region  where  the  absorbed  energy  is  highest.  The  plasma  absorbs  the  incoming  probe 
light,  leading  to  the  dark  area  in  the  center  of  the  melting  region.  This  clearly  emphasizes 
the  need  for  spatial  resolution.  An  integrating  detector  would  give  rise  to  erroneous 
conclusions  about  the  reflectivity.4 

Fig.  3  shows  the  reflectivity  profile  at  the  center  of  the  melting  region,  between 
the  two  white  lines,  in  Fig.  2a.  The  reflectivity  reaches  the  value  for  liquid  silicon  within 
the  30-ps  duration  of  the  pump  pulse.  The  high  time  resolution  of  the  streak  camera 
enables  us  to  confirm  that  the  reflectivity  follows  the  trend  predicted  by  numerical 
simulations  of  heating  above  the  melting  temperature  in  silicon.2  According  to  a  Drude 
model  one  expects  a  decrease  in  reflectivity  of  molten  silicon  when  it  is  heated  above  the 


melting  temperature.  The  laser  fluence  in  Figs.  2a  and  2b  is  470  mJ/cm2  which  is  more 
than  twice  as  large  as  the  melting  threshold  for  silicon  (200  mJ/cm2).  A  numerical 
solution  of  the  one-dimensional  heat  equation  shows  that  at  this  fluence  the  temperature 
of  the  liquid  silicon  exceeds  the  melting  temperature  by  more  than  1000  K. 

Measurements  with  longer  and  more  uniform  probe  pulses  are  currently  in 
progress. 
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1.  Introduction 


Numerous  investigations  on  the  phase  transition  of  silicon  during  picosecond  laser 
annealing  have  been  performed  in  recent  years.  It  has  been  well  established  that  the  silicon 
surface  melts  during  a  picosecond  laser  pulse. 2  Because  liquid  silicon  is  a  metal,  the  reflectivity 
increases  on  melting.  This  has  indeed  been  observed  using  optical  pump-and-probe 
techniques.2-4  Standard  picosecond  pump-and-probe  measurements,  however,  have  some 
serious  inherent  drawbacks.  First,  they  cannot  resolve  reflectivity  changes  that  occur  on  a  time 
scale  of  a  few  picoseconds,  because  they  integrate  over  the  duration  of  the  probe  pulse 
(typically  2(J  ps  or  more).  Second,  they  determine  the  time  profile  of  the  reflectivity  for  every 
pump  fluence  in  a  step-wise  manner  by  varying  the  delay  between  the  pump  and  the  probe 
pulse.  This  introduces  a  large  amount  of  scatter  in  the  data  points,  due  to  shot-to-shot 
variations  in  the  pump  fluence,  and  requires  a  large  amount  of  data  to  be  taken  for  every  time 
profile  of  the  reflectivity.  Also,  they  provide  no  spatial  information  on  the  melting  process. 

To  obtain  spatial  resolution,  a  better  time  resolution,  and  to  measure  the  time  profile  of 
the  reflectivity  on  a  single-shot  basis,  we  use  a  streak  camera  with  a  time  resolution  of  1.8  ps 
for  the  detection  of  the  probe  pulse. 

2.  Experimental  setup 

The  frequency-doubled  output  (30  ps,  532  nm)  of  a  mode-locked  Nd:YAG  laser  is  split 
into  a  probe  and  a  pump  pulse.  The  duration  of  the  probe  pulse  is  stretched  by  splitting  it  in 
four,  delaying  the  four  resulting  pulses  with  respect  to  each  other,  and  recombining  them 
spatially  to  form  a  longer  probe  pulse  of  120  ps  duration.  The  probe  pulse  then  images  an  area 
around  the  100-jim  wide  melting  area  onto  the  entrance  slit  of  a  streak  camera  (Hamamatsu 
Photonics  Cl 587). 

The  entrance  slit  of  the  camera  is  split  into  two  parts.  The  larger  part  is  used  to  image 
the  probe  pulse,  the  smaller  part  to  image  a  fraction  of  the  pump  pulse.  The  latter  acts  as  the 
timing  reference. 

To  enhance  the  sensitivity  of  our  measurements,  the  probe  pulse  is  /^-polarized  and  the 
probe  angle  of  incidence  is  chosen  to  be  65°,  close  to  Brewster’s  angle.  At  this  angle  the 


reflectivity  of  solid  silicon  is  small  (about  10%),  leading  to  an  increase  in  reflectivity  on  melting 
of  a  factor  of  8. 


3.  Results  and  Discussion 

The  streak  images  of  two  measurements  are  shown  in  Figs,  la  and  lb.  The  time  axis  is 
displayed  vertically,  with  time  increasing  from  top  to  bottom.  The  full  height  of  the  image 
corresponds  to  190  ps.  The  horizontal  axis  reflects  the  spatial  profile.  The  pulse  shown  on  the 
left  side,  with  a  width  of  about  30  ps,  is  the  pump  pulse.  The  stretched  probe  pulse  covers 
most  of  the  image.  The  bright  part  is  where  the  silicon  surface  is  melting.  Spatially  it  reflects 
the  Gaussian  intensity  profile  of  the  pump  pulse.  The  time  profile  allows  one  to  study  the 
melting  dynamics. 

Both  measurements  were  performed  at  the  same  laser  fluence.  In  the  case  shown  in 
Fig.  lb,  some  surface  irregularity  leads  to  the  formation  of  a  plasma  on  the  surface  in  the 
region  where  the  absorbed  energy  is  highest.  The  plasma  scatters  the  incoming  probe  light, 
leading  to  the  dark  area  in  the  center  of  the  melting  region.  This  clearly  emphasizes  the  need  for 
spatial  resolution.  An  integrating  detector  would  give  rise  to  erroneous  conclusions  about  the 
reflectivity.5 

Fig.  2  shows  the  reflectivity  profile  at  the  center  of  the  melting  region,  between  the  two 
white  lines,  in  Fig.  la.  The  reflectivity  reaches  the  value  for  liquid  silicon  within  the  30-ps 
duration  of  the  pump  pulse.  The  high  time  resolution  of  the  streak  camera  enables  us  to 
confirm,  that  the  reflectivity  follows  the  trend  predicted  by  numerical  simulations  of  heating 
above  the  melting  temperature  in  silicon.3  According  to  a  Drude  model  one  expects  a  decrease 
in  reflectivity  of  molten  silicon  when  it  is  heated  above  the  melting  temperature.  The  laser 
fluence  in  Figs,  la  and  lb  is  470  mJ/cm2  which  is  more  than  twice  as  large  as  the  melting 
threshold  for  silicon  (200  mJ/cm2).  A  numerical  solution  of  the  one-dimensional  heat  equation 
shows  that  at  this  fluence  the  temperature  of  the  liquid  silicon  exceeds  the  melting  temperature 
by  more  than  1000  K. 

Measurements  with  longer  and  more  uniform  probe  pulses  are  currently  in  progress. 
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ABSTRACT 

We  have  employed  pump-and-probe  techniques  coupled  to 
ellipsometry  to  measure  the  transient  complex  index  of  refraction 
at  1.064  Jim  of  higly  oriented  pyrolitic  graphite  (HOPG) 
interacting  with  20  ps,  .532  Jim  laser  pulses.  When  the  laser  pump 
fluence  exceeds  the  threshold  value  for  melting,  measurement s 
indicate  a  substantial  decrease  of  both  real  and  imaginary  parts 
of  the  index  of  refraction,  thus  confirming  that  molten  graphite 
becomes  less  metallic  Measurements  provide  also  direct  evidence 
of  the  insensitivity  of  our  picosecond  results  to  evaporation 
from  the  irradiated  surface. 


INTRODUCTION 

High  temperature  phases  of  laser  heated  materials  have 
recently  been  studied  with  several  techniques.  In  particular, 
post  irradiation  diagnostics  of  nanosecond  irradiated  samples 
with  optical  and  electron  microscopy,  Raman  scattering, 
Rutherford  backscatter ing  etc.  have  determined  the  existence  of  a 
liquid  state  at  high  temperature  in  layered  compounds  such  as 
graphite  FI, 2].  In  our  laboratory,  instead,  picosecond  time 
resolved  optical  techniques  have  been  used  to  investigate 
directly  Highly  Oriented  Pyrolitic  Graphite  (HOPG)  samples 
during  the  irradiation  with  0.53  Jim,  20  ps  laser  pulses.  The 
advantage  of  this  ultrashort  excitation  is  that  high  surface 
temperatures  are  obtained  with  minimum  irradiation  energy,  while 
optical  probing  at  short  times  avoids  effects  connected  with  the 
evaporation  from  the  sample  surface  [3] . 

Using  time  resolved  pump  and  probe  techniques,  we  have 
previously  reported  a  net  decrease  of  normal  incidence 
reflectivity  for  some  selected  wavelengths  in  a  range  extending 
from  the  IR  to  the  UV,  when  the  laser  excitation  exceeded  a  given 
threshold  fluence  level  [4,5].  The  observations  were  consistent 
with  a  decrease,  at  high  temperature,  of  the  complex  index  of 
refraction  and  thus  of  the  dielectric  function  of  HOPG, 
indicating  a  less  metallic  behavior  of  the  material. 

In  order  to  evaluate  quantitatively  the  transient  complex 
index  of  refraction  nt  =  n  +  i  k  ,  two  simultaneous  measurements 
must  be  performed  on  the  sample  under  identical  conditions.  Since 
HOPG  is  highly  absorbing,  the  easiest  way  is  to  measure  the 
reflectivity  of  the  sample  at  two  different  (orthogonal) 
polarizations  at  various  angles  of  incidence  [6]  .  This  paper 
reports  on  the  application  of  this  ellipsometric  technique  to  the 
case  of  graphite.  The  complex  index  of  refraction  of  HOPG  above 
the  melting  point  has  been  deduced  experimentally  for  a 
wavelength  of  1.06  Jim.  The  qualitative  results  previously 
obtained  [4]  are  confirmed. 
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EXPERIMENTAL  METHOD 

Our  experimental  setup  consisted  of  a  modified  version  of  the 
usual  pump  and  probe  scheme.  HOPG  samples,  with  the  c-axis 
perpendicular  to  the  surface,  were  irradiated  with  0.53  pm,  20  ps 
laser  pulses.  A  second  pulse  at  1.06  pm,  was  passed  through  a 
variable  delay  line  and  a  polarizer  plus  a  X/2  plate,  in  order  to 
control  the  polarization.  The  reflected  beam  was  analysed  by  a 
polarizing  beam  splitter  in  the  two  directions  s-  and  p-  , 
perpendicular  and  parallel  to  the  plane  of  incidence, 
respectively.  Signal  to  noise  ratios  greater  than  10^  for  each 
polarization  channel  were  obtained.  The  sample  could  be  rotated 
so  as  to  achieve  angles  of  incidence  of  the  probe  beam  between 
13°  and  80°.  The  angle  between  pump  and  probe  beams  was  kept 
fixed  at  7°,  a  procedure  which  neither  changed  synchronization 
nor  the  ratio  between  pump  and  probe  spot  areas. 

Absolute  calibration  was  performed  using  glass,  fused  silica 
and  crystalline  silicon  surfaces  as  calibration  standards.  Most 
of  the  errors  associated  with  measurements  illustrated  in  figures 
3-5  are  related  to  spot-to  spot  differences  in  reflectivity  of 
our  HOPG  sample  surface. 


RESULTS  AND  DISCUSSION 

Several  sets  of  experimental  runs  as  a  function  of  incident 
laser  fluence  (J/cm2)  were  performed  with  varying  time  delays 
between  pump  and  probe  pulses  and  at  different  angles  of 
incidence.  The  general  features  of  the  s-ref lectivity  curves 
versus  laser  fluence  are  illustrated  in  figure  1  where  the  probe 
pulse  impinging  on  the  target  at  75°  to  the  normal  is  delayed  by 
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Figure  1:  Graphite  reflectivity  at  1.06  pm  for  s-polarized 
radiation  impinging  at  75°,  as  a  function  of  the  fluence  of  a  20 
ps,  0.53  pm  laser  heating  pulse,  taken  200  ps  after  irradiation. 

200  ps  from  the  pump  pulse.  The  reflectivity  stays  in  general 
almost  constant  at  lower  pump  fluence  levels  and  exhibits  a  slow 
decrease  above  the  critical  fluence  for  melting,  Fth-  This 
decrease  is  more  pronounced  at  low  angles  of  incidence  $  and 
tends  to  level  off  at  high  0  values.  This  behavior  is  similar  to 
what  is  observed  at  normal  incidence.  Absolute  p-ref lectivity 
values  are  always  less  than  the  corresponding  ones  for  the 
s-case.  Their  dependance  versus  laser  fluence  exhibits  a  more 
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complex  behavior.  The  low  fluence  regime  is  characterized  at 
delay  times  just  after  irradiation  by  a  small  reflectivity 
decrease  with  a  recovery  at  fluences  close  to  the  critical  value 
Fth-  This  behavior  seems  consistent  with  the  formation  and  rapid 
decay  of  an  excess  of  plasma  density  near  the  surface,  as  in  the 
case  of  semiconducting  surfaces.  For  d  <  60°,  the  onset  of  the 
phase  transformation  is  indicated  by  a  decrease  of  the 
reflectivity.  For  high  angles  of  incidence,  d  >  60°,  on  the 
contrary,  the  transition  is  marked  by  a  sudden  increase  of  the 
reflectivity,  up  to  a  factor  of  =  2,  which  strikingly  resembles 
what  is  observed  in  the  case  of  semiconductors  like  silicon. 
Figure  2  illustrates  one  such  case  for  0  =  75°  and  a  delay  of  200 
ps .  As  for  silicon  a  reduction  of  the  reflectivity  is  observed  at 
high  fluence.  This  may  be  interpreted  as  due  to  overheating  of 
the  absorbing  liquid  surface  and  to  the  effects  of  surface 
evaporation  and  ionization  in  a  similar  manner  as  occurred  in 
liquid  silicon. 
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Figure  2:  Graphite  reflectivity  at  1.06  Jim  for  p-polarized 
radiation  impinging  at  75°,  as  a  function  of  the  fluence  of  a  20 
ps,  0.53  Jim  laser  heating  pulse,  taken  200  ps  after  irradiation. 

Other  information  on  these  processes  may  be  obtained  by 
plotting  our  reflectivity  data  as  a  function  of  time  for 
different  excitation  levels.  In  particular,  figure  3  shows  the 
time  evolution  of  s-ref lectivity  at  75°  incidence  for  fluences  at 
and  above  the  critical  valve  Fth  for  surface  melting.  For  all 
three  fluence  levels  illustrated,  the  reflectivity  appears  to 
decrease  initially  during  the  laser  pulse.  At  F  =  Fth,  the 
decrease  is  followed  by  a  recovery  of  the  reflectivity  in  =  200 
ps  to  a  constant  final  value  Ru.  At  higher  fluences  the 
intermediate  reflectivity  variations  disappear.  Instead,  a 
sharper  drop,  continuing  up  to  times  =  500  ps  after  irradiation, 
is  apparent.  It  is  followed  by  a  slow  recovery  to  R<»,  presumably 
connected  with  the  cooling  of  the  surface. 

The  corresponding  variations  of  p-ref lectivity  with  time  are 
shown  in  figure  4.  Here,  the  increase  in  reflectivity  occurring 
at  near-threshold  fluence  values  only  at  later  times  indicates 
that  for  the  first  200  ps  after  irradiation  the  cloud  of 
evaporating  material  is  not  responsible  for  the  complex 
reflectivity  variations  observed,  which  we  attribute  to  the  hot 
HOPG  surface.  The  characteristic  times  for  reflectivity  increase 
are  clearly  dependent  upon  laser  input  fluence.  They  are 
consistent  with  an  ultrafast  phase  transformation  occurring  at 
earlier  times  with  higher  irradiation  levels.  The  fact  that 
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Figure  3:  Time  evolution  of  s-ref lectivity  in  graphite  for  1.06 
Jim  radiation  at  75°  incidence  angle,  for  the  20  ps,  0.53  Jim  pump 
fluence  levels  indicated.  The  critical  fluence  for  surface 
melting,  Fth,  corresponds  to  0.14  J/cm2. 


Figure  4:  Time  evolution  of  p-ref lectivity  in  graphite  for  1.06 
Jim  radiation  at  75°  incidence  angle,  for  the  20  ps,  0.53  Jim  pump 
fluence  levels  indicated.  Note  the  different  scale  on  the 
reflectivity  axis. 

the  maximum  of  reflectivity  is  reached  long  after  the  exciting 
laser  pulse  may  be  due  to  the  presence  of  interference  effects 
connected  with  the  formation  and  evolution  of  a  molten  layer  on 
the  surface  of  the  sample.  The  observed  changes  in  p-ref lectivity 
reveal  the  extremely  fast  dynamics  of  the  heating  and  cooling  of 
HOPG.  The  decay  of  the  high  reflectivity  phase  can  be  interpreted 
as  the  onset  of  surface  cooling  due  to  both  thermal  conduction 
into  the  bulk  and  evaporation  of  material  from  the  surface. 

The  variation  of  the  experimental  data  taken  at  different 
angles  of  incidence  are  summarized  in  figure  5,  where 
experimental  p-  and  s-ref lect ivity  obtained  with  F  =  1.75  Fth, 
=200  ps  after  HOPG  irradiation  are  plotted  versus  angle  of 
incidence.  The  continuous,  thin  lines  in  figure  5  denote 
reflectivity  of  unirradiated  HOPG  as  derived  from  our  calibration 
measurements.  These  results  correspond  to  an  initial  ordinary 
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complex  index  of  refraction  nc  =  2.9  +  i  2.06,  in  agreement  with 
published  cw  data  [7]  .  Since  graphite  is  a  negative  uniaxial 
crystal,  its  extraordinary  index  has  some  influence  on  the 
P~ re f 1 ec t i v i t y  of  the  HOPG  sample.  The  low  temperature 
p-ref lectivity  can  be  fitted  best  by  assuming  ne  =  2  .  +i  0.25 
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Figure  5:  1.06  pm  s  (upper  curves)  and  p  (lower  curves) 
reflectivities  of  graphite  versus  angle  of  incidence.  Continuous 
thin  line:  best  fitted  values  at  room  temperature,  evaluated  with 
n0  =  2.9  +  i  2.06,  .  Data  points:  measured  reflectivities  =200  ps 
after  a  20  ps,  0.53  pm,  1.75  Ft excitation  laser  pulse.  The 
thick  continuous  curve  is  calculated  from  the  best  fitting  index 
of  refraction  nht  =  1.3  +  i  1.2. 


The  two  curves  for  p-ref lectivity  versus  angle  of  incidence 
for  room  and  high  temperature,  respectively,  cross  each  other  at 
an  angle  between  60  and  70°.  This  implies  that  the  pseudo 
Brewster  angle  i3HT  for  the  high  temperature  case  is  smaller  than 
the  corresponding  angle  dLT  at  room  temperature.  This  is  per  se 
an  indication  that  graphite  has  a  lower,  transient  index  of 
refraction  when  excited  above  the  critical  fluence  level  Fth. 
The  reflectivity  measurements  performed  in  this  time  range 
obviously  are  ilqlL  affected  by  optical  opacity  generated  by  an 
evaporating  cloud  of  material,  since  reflectivity  changes  at  high 
fluence  have  an  opposite  sign  with  polarization. 

Assuming  isotropic  conditions  for  the  high  temperature 
medium,  the  data  illustated  in  figure  5  are  consistent  with  a 
complex  index  of  refraction  nht  =  1.3  +  i  1.2,  evaluated  by  best 
fitting  to  the  experimental  points.  This  value  for  n  is  smaller 
than  the  initial  one,  in  agreement  with  our  earlier  measurements 
and  interpretation  [4].  In  particular  the  pseudo  Brewster  angle 
in  the  high  temperature  phase  is  smaller  than  for  room 
temperature  HOPG,  as  the  increase  of  p-ref lect ivity  for  high 
angles  of  incidence  indicates.  The  data  indicate  that  the  high 
temperature  phase  is  less  metallic  in  character  than  graphite, 
and  may  even  exhibit  a  band  gap.  The  layered  structure 
characteristic  of  graphite  disappears  due  to  cross  linking. 
Perhaps  the  high  temperature  phase  could  be  better  modeled  as  a 
distorted  diamond  structure.  Theoretical  band  calculations,  based 
on  an  isotropic  (cubic)  structure,  predict  the  existence  of  a 
bandgap  (8) .  Systematic  investigations  over  a  wavelength  interval 
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extending  from  the  ultraviolet  to  the  infrared  may  eventually 
reveal  the  electronic  structure  of  the  high  temperature  phase. 
Such  experiments  are  feasible  by  the  techniques  described  in  this 
note  . 


CONCLUSIONS 

By  coupling  ultrashort  laser  techniques  with  ellipsometry, 
the  index  of  refraction  of  materials  in  transient,  out  of 
equilibrium  conditions  may  be  directly  evaluated.  In  particular, 
time  resolved  optical  properties  of  materials  heated  above 
melting  may  be  studied  before  substantial  evaporation  occurs.  By 
extending  this  technique  to  a  wide  wavelengths  interval, crucial 
information  on  the  electronic  structure  of  molten  graphite  might 
be  obtained. 


ACKNOWLEDGEMENTS 

We  wish  to  thank  Dr  .  Chao  Yen  Huang  for  many  important 
comments  and  suggestions  throughout  this  work,  which  has  been 
supported  by  the  Office  of  Naval  Research  under  contract 
N.00014-84-K-04654 


REFERENCES 

1.  T .Venkatesan,  D . C . Jacobson,  J.M. Gibson,  E.S. Elman,  G.Braunstein 
M . S . Dresselhaus  and  G.  Dresselhaus,  Phys . Rev . Letters  51,  360 
(1984) 

2.  J . S . Speck, J . Steinbeck,  G.Brunstein,  M. S . Dresselhaus,  and 
T. Venkatesan,  Mat.  Res.  Soc.  Symp.Proc.  51,  263  (1986) 

3.  N .  Bloembergen,  Mat.  Res.  Soc.  Symp.Proc.  51,  3  (1986) 

4.  A.M.Malvezzi,  N . Bloembergen  and  C.Y. Huang,  Phys . Rev . Letters  57 . 
146  (1986) 

5.  C.Y. Huang,  A.M.Malvezzi,  N . Bloembergen  and  J.M. Liu,  Mat.  Res. 
Soc.  Proc.  51,  245  (1986) 

6.  P.M.Fauchet,  W.L.Nigham  Jr,  Appl .  Phys.  Letters  IS,  785  (1986) 

7.  E. A. Taft  and  H.R. Philipp,  Phys . Rev . Hfi,  A197  (1965) 

8.  M. Cohen,  private  communication 


269 


TIME-RESOLVED  PICOSECOND  REFLECTIVITY  STUDY  OF  LASER-EXCITED 

LAYERED  COMPOUNDS 


C.Y. HUANG*,  A . M . MALVEZZ I  *  * ,  N.  BLOEMBERGEN *  *  AND  F.J.DI  SALVO, 
*  *  * 

JR. 

Los  Alamos  National  Laboratory,  Los  Alamos,  New  Mexico  USA 
Division  of  Applied  Sciences,  Harvard  University,  Cambridge, 
MA  02138,  USA 

***  AT&T  Bell  Laboratories,  Murray  Hill,  N.J.  01  974,  USA 


ABSTRACT 

We  have  employed  the  pump-and-pr obe  technique  to  perform 
picosecond  time  resolved  measurement  s  of  the  reflectivity 
changes  in  two  archetypal  layered  compounds,  lT-TiS2  and 
lT-TiSe2  probed  at  1.064  Jim  after  pumping  by  20  ps,  .532  Jim 
laser  pulses.  At  the  threshold  fluence,  -  40  mJ/cm2,  the 
reflectivity  drops  sharply,  marking  the  occurrence  of  a  phase 
transformation  on  the  surface  of  the  sample.  Above  threshold, 
the  reflectivity  reaches  a  value  as  low  as  -0.1  at  high 
fluences,  strongly  suggesting  that,  like  in  graphite,  the  high 
temperature  phase  is  not  metallic. 


INTRODUCTION 

Very  recently,  the  transient  optical  properties  of 
crystalline  graphite  surfaces  excited  by  visible  picosecond 
laser  pulses  have  been  investigated  with  20  -  30  ps  time 
resolution  1#2.  Reflectivity  measurements  at  low  excitation 
levels  have  revealed  the  presence  of  a  fast  decaying 
elecron-hole  plasma  evolving  from  the  initial  semimetallic 
configuration.  More  interestingly,  at  a  well  defined  pump 
threshold  value  of  -140  mJ/cm2  a  phase  transition  occurs,  which 
is  characterized  by  lower  reflectivity  values.  We  have  concluded 
that  the  high  temperature  phase  is  not  metallic,  having  a  lower 
value  of  the  index  of  refraction,  but,  rather,  intermediate 
between  the  one  of  graphite  and  that  of  diamond.  This  unexpected 
new  phase  is  presumably  caused  by  a  strong  perturbation  of  the 
spatial  arrangement  of  carbon  atoms  in  the  graphite  lattice, 
which  consists  of  hexagonally  arranged,  loosely  stacked  atomic 
layers.  The  question  thus  arises  whether  the  electronic 
properties  of  the  high  temperature  phase  originate  from  a 
particular  initial  structure.  For  this  reason,  it  is  of  interest 
to  investigate  other  archetypal  layered  compounds. 

In  this  paper,  we  have  chosen  two  crystalline  layered 
compounds,  lT-TiS2  and  lT-TiSe2 .  The  titanium  atoms  are 
octahedrally  coordinated  with  chalcogen  atoms  forming  a 
chalcogen-metal-chalcogen  sandwich  layer  configuration.3  In  a 
plane,  atoms  are  hexagonally  packed.  However,  the  coordination 
around  the  non-metal  is  quite  lopsided,  leading  to  the  marked 
cleavage  properties  perpendicular  to  the  hexagonal  symmetry 
axis.  The  basic  atomic  structure  of  loosely  coupled  X  -  Ti  -  X 
(X  =  S,  Se)  atom  sheet  sandwiches  makes  the  mechanical,  thermal 
and  electrical  properties  extremely  anisotropic,  like  those  of 
graphite.  The  lattice  constants  for  TiS2  4  are  slightly  smaller 
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than  these  of  TiSej.5  It  is  now  well  established  that  the  former 
is  an  extrinsic  semiconductor  with  an  indirect  band  gap  of  0.2  - 
0.3  eV,  6-9  but  it  is  normally  non-stoichiomet  r  ic  and  thus 
exhibits  metallic  conductivity  arising  from  an  excess  of  Ti 
atoms.  TiSeo  is  an  intrinsic  semimetal.  0-12  All  these  physical 
properties  are  similar  to  those  of  graphite.  Hence,  it  is  of 
importance  to  perform  experiments  on  the  high  temperature  regime 
for  these  compounds  and  then  compare  the  results  with  those  of 
highly  oriented  pyrolitic  graphite  (HOPG).1'2 


EXPERIMENTAL  RESULTS 

We  have  performed  a  series  of  time-resolved  picosecon 
pump-and-probe  measurements  similar  to  those  made  on  HOPG.  1( 
The  pumping  was  provided  by  20  ps  pulses  at  0.532  |im  by 
frequency  doubling  the  output  of  an  active-passive  mode-locked 
Nd-YAG  laser  system.  The  variable  time  delayed  1.064  pm  pulses 
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were  used  to  monitor  the  reflectivity  at  the  center  of  the 
excited  surface.  Near  normal  incidence  angles  are  employed  for 
both  pump  and  probe  laser  beams.  All  the  samples  used  in  this 
work  were  single  crystals. 

We  have  found  that  a  damage  spot  appears  on  the  sample  above  a 
fixed  green  pump  fluence  value.  Our  optical  micrographs  show  the 
appearance  of  a  black  ring  marking  the  damage  area.  Its 
formation  suggests  the  occurrence  of  a  phase  transformation.  The 
size  of  the  ring  strictly  follows  the  spatial  variations  of  the 
gaussian  fluence  distribution  of  the  laser  beam,  indicating  the 
thermal  origin  of  the  effect.  The  extrapolated  value  of  the 
threshold  fluence  is  Fth  =  43  mJ/cm2 

Figure  1  shows  our  reflectivity  results  obtained  on  lT-TiS2 
single  crysLal  probed  at  1.064  pm  as  a  function  of  the  0.532  pm 
pumping  fluence  at  the  different  delay  times  td.  As  Figure  1  (a 
-  d)  shows,  the  reflectivity  drops  sharply  around  40  mJ/cm2  fcr 
all  positive  delay  times  td,  in  agreement  with  the  value  of  Fth 
obtained  from  the  post-experimental  examination  of  the  damaged 
-,rea  discussed  above.  Consequently,  like  in  HOPG,  we  may 
identify  Fth  as  the  threshold  fluence  for  an  ultrafast  phase 
transition.  The  drop  is  particularly  pronounced  at  =100  ps, 
where  the  reflectivity  reaches  a  value  as  low  as  =  0.1  from  the 
initial  value  of  0.52.  This  variation  exceeds  the  similar 
variation  in  the  reflectivity  observed  in  HOPG  at  the  same 
wavelength.  It  has  been  shown  that  for  a  delay  time  of  =  10  ps, 
only  one  atomic  layer  can  escape  the  surface  of  the  sample,13 
giving  a  negligible  contribution  to  the  reflectivity.  As  shown 
in  Fig.  l-(a),  even  when  the  probe  pulse  overlaps  the  pump  pulse 
(td  -  0  ps)  the  reflectivity  drops  above  the  threshold  fluence 
value,  thus  clearly  demonstrating  that  the  decrease  of  the 
reflectivity  is  not  due  to  evaporation. 


Figure  2.  Reflectivity  of 
1  T - T i S 2  at  1.064  ^m  as  a 
function  of  probe  delay 
time  td,  for  various  pump 
fluences.  The  origin  of  td 
is  shifted  by  50  ps . 


In  Figure  2  we  plot  the  measured  reflectivity  as  a  function 
of  the  delay  time,  td.  Note  that  the  time  origin  in  the  abscissa 
has  been  shifted  by  50  ps  in  order  to  accommodate  the  data  taken 
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before  the  arrival  of  the  pump  pulse.  At  25  mJ/cm2,  the 
reflectivity  is  practically  time  independent.  At  fluences 
greater  than  Fth,  the  reflectivity  values  are  decreasing  at 
first,  reaching  a  minimum  at  =  200  ps,  and  then  recover  to 
higher  levels.  This  behavior  of  the  reflectivity  clearly 
indicates  the  occurrence  of  a  transient  new  phase  induced  by  the 
pump  laser  pulse  and  presumably  terminated  by  heat  dissipation 
in  the  bulk  of  the  crystal.  The  lack  of  data  on  the  thermal 
properties  of  these  compounds  prevents  quantitative  estimates. 
The  final  value  =  .35  of  the  reflectivity  taken  at  very  long 
delay  times  is  lower  than  the  initial  one,  indicating  the 
presence  of  a  permanent  damage  on  the  surface  as  soon  as  the 
threshold  fluence  value  Ft^  is  exceeded. 


Figure  3.  Reflectivity  of 
lT-TiSe2  at  1.064  |im  as  a 
function  of  0.532  urn  pump 
fluence  at  the  delay  (a)  - 
0  ps  and  (b)  -  120  ps. 


0  100  200  300  400  300 
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The  reflectivity  data  for  the  semimetal  lT-TiSe2  , 
illustrated  in  Figure  3,  are  qualitatively  similar  to  those  for 
the  semiconductor  lT-TiS2.  Also,  the  temporal  behavior  at  the 
different  fluences  appears  similar  to  the  one  presented  in  Fig 
2.  Our  results  indicate  that  the.  high  temperature  optical 
behavior  of  these  two  materials  is  independent  of  the  initial 
electronic  properties  of  these  two  layered  compounds. 


CONCLUSIONS 

The  transient  optical  properties  of  TiS2  and  TiSe2  crystals 
under  picosecond  visible  excitation  have  been  investigated  at 
1.064  ^m  with  picosecond  time  resolution.  Below  the  threshold 
fluence  of  =40  mJ/cm2  for  20  ps,  0.532  urn  pump  laser  pulses,  no 
significant  changes  in  the  reflectivity  at  1.064  Jim  for  both 
materials  have  been  observed.  Due  to  the  small  value  of  the  band 
gap  in  TiS2  as  compared  to  the  probing  photon  energy,  the  plasma 
contribution  to  the  dielectric  constant  appears  to  be 
negligible.  Presumably,  longer  wavelengths  and/or  larger 
incidence  angles  have  to  be  used  to  detect  plasma  effects  from 
the  optical  response  of  these  materials. 
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At  fluences  exceeding  the  threshold  value  Fth  the 
reflectivity  decreases  to  values  as  low  as  =  0.1.  These 
reflectivity  changes  occur  during  the  exciting  laser  pulse.  This 
suggests  that  the  high  temperature  phase  of  these  compounds  is 
not  metallic.  The  laser  pump  pulse  gives  rise  to  a  transition 
from  a  semimetallic  or  small-band-gap  semiconducting  electronic 
structure  into  a  rather  transparent  phase,  where  a  smaller 
number  of  carriers  interacts  with  the  incoming  laser  pulse. 
Light  can  therefore  travel  into  the  material  for  distances 
exceeding  the  original  extinction  depth.  Comparison  14  of  the 
morphology  of  picosecond  with  nanosecond  laser  irradiated 
graphite  indicate  that  ultrafast  melting  occurs  at  fluences 
exceeding  Fth .  Similar  considerations  apply  here  for  TiSe2  and 
TiS2  ,  the  critical  value  for  melting  being  lower  due  to  the 
different  optical  and  thermal  properties  of  these  compounds. 

This  comparative  study  suggests  that  a  similar,  non-metallic 
high  temperature  phase  exists  for  archetypal  layered  structures 
other  than  graphite.  The  electronic  properties  at  high 
temperatures  could  originate  from  the  initial  presence  of  a 
layered  structure  which  undergoes  heavy  disruption  under  laser 
excitation.  A  deeper  insight,  however,  is  necessary  to  unravel 
the  true  nature  of  the  structure.  The  transient  optical 
properties,  as  derived  i.e.  by  pump-and-probe  experiments  or 
time-resolved  ellipsometry  at  different  wavelengths,  can  provide 
useful  information  if  coupled  with  other  firm  optical  and 
thermal  data  on  the  materials  under  study. 
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1  Int-OdUCt'O^ 

The  irradiation  of  a  strongly  absorbing  solid  surface  with  a  short  laser 
ouise  is  a  means  of  concentrating  thermal  energy  in  a  very  thin  layer.  The 
timescale  considered  in  this  paDer  is  greater  than  1  ps,  which  is  longer 
than  the  time  constants  for  most  electronic  relaxation  processes,  so  that 
only  thermal  effects  need  be  considered  [1].  In  a  strongly  aDsoroing  solid 
(metal,  UV  on  semiconductors)  the  absorDtion  depth,  a*1  (a:  aosorption 
coefficient),  'SJess  than  the  thermal  diffusion  length  during  irradiation, 
!t  =  (D.htp)'/2  (D»h:  thermal  diffusivitv;  tp :  pulse  length),  F or  a  30  ps 
Dulse  oH  [ron,  the  'example  used  in  all  further  discussion,  these  lengtns 
are  200  A  and  500  A,  respectively.  A  general  discussion  of  the  two 
regimes  corresponding  to  the  relative  values  of  a"'  and  ly  can  be  found  in 
a  reytew  paoer  by  Bloembergen  [2].  It  is  worth  noting  that  for  t0  <  Ips.  ly 
<  a‘ 1 ,  so  that  shortening  the  pulse  further  no  longer  leads  to  a  steeper 
temperature  gradient  and  a  higher  quenc.n  rate. 

If  the  fluence  is  low  enough  so  that  r\p  visible  permanent  surface 
damage  occurs  (typically  around  100  mJ/cm2),  a  liquid  layer  of  thicKness 
ly  with  a  temperature  of  several  thousand  degrees  is  created  during  the 
irradiation.  Further  melt-in,  to  a  maximum  thicKness  d  of  a  few  thousand 
A.  occurs  until  the  crystal-melt  interface  temperature  drops  to  the 
eauiiibnum  melting  temperature,  Tm;  this  stage  lasts  on  the  order  of  ICO 
ds.  Upon  further  cooling,  the  interface  temperature  drops  below  Tm,  and 
the  crystal  regrows.  The  thermal  parameters  of  this  cooling  process  can 
be  estimated  simply  by  a  dimensional  analysis  based  on  the  thicKness  of 
the  molten  layer,  d,  and  the  temperature  of  the  melt,  Tm,  as  snown  in 
Taple  1  [3]. 


Table  1.  Thermal  Parameters  in  Melt  Quenching  of  Metals 


with  a  30  ps  Laser  Pulse. 

From  [3] 

Melt  temperature 

Tm 

-  IQ3  K 

Melt  thickness 

d 

-  10‘7  m 

Temperature  gradient 

VT  -  Tpp/d 

-  1010  K/m 

Cooling  rate 

t  -  Dth  VT/d 

-  1012  K/  S 

Melt  lifetime 

x  m  Ty/T 

-  10’9  s 

Isotherm  velocity 

Uy  -  t/VT 

*  100  m/s 

Heat-flow  limited 
crystal  growth  velocity 

U(-|  ■  k  VT/  AHc  v 

-  230  m/s 
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The  metastable  phases  formed  in  this  process  are:  the  crystal 
overheated  above  T_,,  the  liquid  overheated  above  the  boiling  point,  Tw, 
the  liquid  undercooled  below  Tm ,  and,  possibly,  new  metastabie 
crystalline  phases  or  glasses.  This  paper  gives  a  concise  review  of  the 
thermodynamic  stability  of  these  phases,  and  the  kinetics  of  the 
transitions  between  them. 

2  Kine>;cs  of  Crystallization 

The  theory  of  the  motion  of  interphase  boundaries  has  been  reviewed  in 
detail  m  a  number  of  papers  [4],  The  results  for  crystallization,  which  are 
easily  generalized  for  other  transitions,  can  be  summarized  as  follows. 
The  crystal  growth  velocity  can  be  written  as 

u  -  f  k  X  [1  -  exp(-  ]  ,  (1) 

where  k  is  the  atom  jump  frequency  across  the  crystal-melt  interface,  f 
is  the  fraction  of  interface  sites  that  can  incorporate  a  new  atom,  X  is 
the  interatomic  distance,  aGc  is  the  difference  in  molar  free  energy 
between  crystal  and  melt,  and  Tj  is  the  interface  temperature.  If  the 
undercooling  of  the  interface,  AT,  -  Tm-Tj,  is  small,  (1)  can  be  linearized 
to 


u  =  k  X  (  at,/  Tj)  .  (2) 

The  crystal  growth  rate  also  depends  on  the  rate  at  which  the  heat  of 
crystallization,  aHc  „  per  unit  volume,  can  be  removed.  Since  this  must 
occur  by  heat  flow 'down  the  temperature  gradient  at  the  interface,  V , T , 
the  neat  flux,  J,  must  ooey  the  relations 

J  =  u  aHc  v  =  k  VjT  ,  (3) 

wnere  vis  the  thermal  conductivity.  Equation  r 5 j  can  be  turned  around  to 
define  a  heat  flow-limited  velocity: 

U(-,  =  k  VjT/  aHc  v  ,  (4) 

whicn.  for  the  example  used  here,  is  also  listed  in  Table  1.  Combining  (2) 
anq  (4)  then  gives  for  metais 

Uh./  ^  ^  x  (Tm  *  Tj)/  Tj  •  (5) 

The  maximum  growth  velocity  of  a  crystal,  umaX-  k.X ;  is  determmea 
by  the  nature  of  the  atomic  rearrangements  that  control  the  growtn 
process.  In  pure  metais,  dilute  alleys  or  comoounds  with  simple  crystal 
sturctures.  k  can  be  taken  as  the  frequency  of  thermal  vibration  (-  Deove 
‘recuencyl,  so  that  ur-Lax  -  us,  the  speed  of  sound  in  the  liquid,  which,  'is 
STrcu.nc  3000  m/s  for  "re.  This  is  called  the  collision-controlled  regime, 
because  the  growth  occurs  by  simple  im.omgement  of  the  liquid  atoms 
onto  the  crystal.  This  is  confirmed  by  direct  measurements  of  the  growth 
velocity  of  pure  metais  and  disordered  alloys  (Au,  Cu,  Cu-Au)  by 
transient  reflectivity  [5]  which  yield  large  values  (>100  m/s)  that  car, 
only  be  accounted  for  by  a  collision-controlled  mechanism.  From  (5)  and 
the  values  from  Table  I,  the  interface  in  Fe  is  found  to  be  undercocled  by 
only  128  K.  An  undercooling  of  at  least  295  K  is  required  for  homogeneous 
crystal  nucleation  to  be  appreciable  [6]. 

In  more  concentrated  alloys  or  in  compounds  with  a  more  complicated 
crystal  structure,  where  a  change  of  the  nearest  neighbor  configuration, 
i.e.  diffusion,  is  necessary,  the  jump  frequency  can  be  taken  as  D/a.2,  and 
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umax  D  =  Di/X,  which  is  on  the  order  of  10  m/s.  This  is  the 
ditfusFon-controiled  regime.  Table  I  shows  that  for  picosecond  laser 
quenching  Umax  q  <<  uj,  up,  so  that  this  type  of  growth  can  easily  be 
suppressed.  Tne  Interface  temperature  drops  far  below  Tm,  and  glasses 
can  be  formed.  This  is  confirmed  by  numerous  experiments  on  binary 
metallic  systems  (eg.,  Ni-Nb  [7]),  in  which  glasses  have  been  formed  by 
ps  irradiation  at  compositions  where  tne  intermetallic  crystalline 
compounds  have  a  large  unit  cell,  or  in  which  long-range  chemical  order 
must  be  established. 


Melting  of  a  crystal  is  a  first-order  transition,  which  usually  occurs 
heterogeneously,  i.e.,  by  nucieation  of  the  liquid  and  motion  of  a 
liquid-crystal  interface.  Since  for  most  materials  ocv  >  0|v  +  Oq|  (oji: 
surface  tension  of  the  interface  between  phases  i  and  j;  c:  crystal;  v: 
vapor;  I:  liquid),  there  is  no  barrier  to  nucieation  of  the  liquid  at  the 
surface  of  the  pure  crystal  [8].  It  should  be  kept  in  mind,  howev--r,  that,  if 
the  crystal  is  covered  with  a  high  melting  layer  (e.g.,  an  oxide)  and  the 
surface  tension  of  the  resulting  interface  is  low,  nucieation  of  the  melt 
may  require  some  overheating  of  the  crystal. 


The  motion  of  the  crystal-melt  interface  is  governed  by  the  kinetics 
described  in  Section  2.  During  the  initial  stages  of  the  melt-in  phase  the 
temperature  gradients  are  an  order  of  magnitude  greater  than  during 
regrowth.  This  means  that  the  heat-flow  limited  crystal  growth  velocity, 
given  by  (4)  is  also  an  order  of  magnitude  greater  than  the  value  listed  in 
Table  I,  or  about  2000  m/s.  For  pure  metals,  where  Upp,^  =  us  ~  3000 
m/s,  (5)  gives  for  the  overheating:  ATj=  0.66  Tj,  which  for  Fe  corresponds 
to  an  overheating  of  several  thousand  degrees.  For  materials  where 
melting  is  diffusion-limited,  the  overheating  is  even  higher. 


The  overheating  required  to  make  homogeneous  nucieation  of  the  melt 
(i.e.,  in  the  interior  of  the  crystal)  the  dominant  mechanism  can  be 
estimated  from  classical  nucieation  theory  [8],  which  gives  the  following 
expression  for  the  steady  state  nucieation  frequency  (number  of  nuclei 
per  unit  volume  and  time): 


I  -  N  k  exp  (  - 


(6) 


where  N  is  the  number  of  atoms  per  unit  volume,  k  is  the  interfacial  jump 
frequency,  aGv  is  the  free  energy  difference  per  unit  volume,  o  is  the 
crystal-melt  interfacial  tension,  and  kg  is  Boltzmann's  constant. 


The  minimum  nucieation  frequency  requirea  can  be  estimated  as:  Ipp j n 
=  a  /  a2  tD,  where  a  is  the  minimum  lateral  distance  between  nuclei  for 
inducing  routing  throughout  the  layer.  Using  a  «  1pm  gives  a  value  for 
ImjP,  of  1030  m'Js'T.  (It  should  be  kept  in  mind  that  the  analysis  is  very 
insensitive  to  the  choice  of  lmjn).  Using  a  conservatively  high  estimate 
for  c,  0.3  N/m  [6j,  an  overheating  of  600  K  is  found  to  be  sufficient  to 
exceed  lmjn  in  he.  The  transient  time  for  establishing  the  steady-state 
nucieation  fate  [9]  can  be  estimated  as  the  time  required  to  assemble  a 
critical  nucleus  of  radius  r*»  -2o/aGv,  which  is  r*/kX;  for  the  example  of 
Fe  overheated  by  600  K,  this  is  only  .a  few  ps,  so  that  homogeneous 
nucieation  of  the  melt  can  be  considered  possible  here. 


For  systems  where  the  melting  kinetics  are  relatively  slow  (i.e.  the 
jump  frequency,  k,  scaling  both  nucieation  and  growth  is  small)  it  is 
possible  to  overheat  far  enough  to  reach  the  stability  limit  of  the  crystal. 
This  is  the  temperature  at  which  the  shear  modulus,  p,  vanishes,  and  at 
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Fig.  1 .  P-T  diagram  for  Fe. 
calculated  from  the  van  cer 
Waals  model  and  the 
critical  parameters 


wr.icn  melting  can  occur  homogeneously  by  unlimited  thermal  shear 
emplacements  throughout  the  crystal.  An  extrapolation  based  on  typical 
va.ues  of  op'dT  [10]  sets  this  staoility  limit  at  about  2Tm. 

4  The  Qve^ea ted  Jquid 


That  licuics  can  be  overheated  considerably  above  their  melting 
temoerature  is  a  well-known  experimental  fact  [11],  which  has  been 
Qualitatively  understood  since  van  der  Waais  (Figure  1).  The  liquid  can  be 
overneatec  into  the  metastaole  region,  wnere  nucleation  of  a  vapor  bubble 
is  necessary  to  induce  the  transformation,  up  to  the  spmodal  line,  where 
tne  liquid-vapor  surface  tension,  C|V,  and  hence,  the  nucleation  barrier, 
vanish.  Figure  1  shows  that  at  atmospheric  pressure  (0.1  MPa),  this 
instability  occurs  at  5700  K. 

A  calculation  of  the  rate  of  homogeneous  bubble  nucleation.  based  on 
the  classical  expression  similar  to  (6),  ta*mg  into  account  the  roughiy 
unear  cecrease  of  civ  wun  temperature  to  the  spinocal  [1 2] ,  shows  that 
the  temperature  of  the  overheated  liquid  must  be  within  a  few  huncreo 
degrees  from  the  spmodal  for  the  nucleation  rate  to  exceed  l^,,n  =  10°^ 
m  '^s  ,  .estimated  above.  An  estimate  of  the  transient  time  for  this  . 
transformation  is  r*/r  X **  (r  -  p/VZurnkT,  the  impingement  rate  of  the 
vapor  atoms),  which  is  7  ps  for  the  present  example.  Given  that  this  is  a 
lower  limit  (perfect  sticking  is  assumed  in  calculating  D,  and  that  the 
extreme  overheat  required  Tor  nucleation  is  only  present  during  a  fraction 
of  tD,  the  transient  effects  probably  enhance  the  stability  of  the 
ovemeated  liquid. 

In  a  typical  ps  melt  quenching  experiment  the  fluence  is  not  high 
enough  to  reach  the  extreme  temperatures  required  for  bulk  vaporization. 
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Only  surfacs  evaporation  needs  to  be  consicered,  which  can  be  estimated 
conveniently  from  the  impingement  rate,  r,  of  the  equilibrium  vapor 
phase.  It  has  been  shown  that  the  loss  of  mass  and  energy  by  this  mode  is 
negligible  in  ps  experiments  [13]. 

The  permanent  surface  damage  (holes)  sometimes  observed  in  ps 
experiments  is  therefore  not  due  to  direct  evaporation,  but  to  mechanical 
displacement  of  the  liquid  by  the  recoil  pressure  of  the  evaporating 
atoms.  This  is  confirmed  by  the  morphology  of  the  displaced  material. 
Direct  evidence  of  the  recoil  pressure  can  be  obtained  by  observing 
the  vibration  of  a  metallic  reed  pulse  irradiated  at  one  end  and  clamped  at 
the  other  [14], 

The  author's  work  in  this  area  is  supported  by  the  Office  of  Naval  Research  under 
contract  number  N0001 4-85-K-0634. 
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METASTABLE  STATES  IN  PULSED  LASER  QUENCHING 


Frans  Spaepen 

Division  of  Applied  Sciences 
Harvard  University 
Cambridge,  MA  02138 


It  Is  reviewed  how  irradiation  of  a  metallic  surface  with  a  short  laser  pulse 
leads  to  the  successive  formation  of  these  metastable  states:  overheated 
crystal,  overheated  liquid,  undercooled  liquid,  and  metastable  crystal  or 
glass.  The  stability  of  a  liquid  overheated  above  Its  melting  temperature  is 
analyzed.  The  thermodynamics  and  kinetics  of  melting  and  crystallization  of 
pure  metals  and  alloys  are  reviewed.  The  results  of  quenching  experiments  on 
the  Nb-Si  system  are  used  to  illustrate  how  phases  with  high  isotropy  and 
small  unit  cell  size  are  kinetically  favored. 


Introduction 


Irradiating  the  surface  of  a  strongly  absorbing  solid  (metal,  semicon¬ 
ductor)  with  a  short  laser  pulse  is  a  method  for  depositing  thermal  energy  in 
a  very  thin  surface  layer,  which,  as  a  result,  is  first  melted  and  subsequently 
quenched  at  a  very  high  rate.  Picosecond  laser  pulsing  can  produce  quench 
rates  exceeding  10^  K/s,  making  It  the  fastest  liquid  quenching  method 
available. 

In  the  first  part  of  this  paper  the  three  stages  of  the  process  (energy 
deposition,  melt-in,  and  cooling)  are  analyzed  in  some  detail,  and  three  meta¬ 
stable  states  are  identified:  the  overheated  crystalline  solid  (above  its 
melting  point),  the  overheated  1 iquid  (above  the  boiling  point),  and  the 
undercooled  liquid  (below  the  melting  point). 

The  second  part  of  the  paper  is  a  detailed  analysis  of  the  overheated 
liquid  state.  It  is  shown  that  the  liquid  can  remain  in  this  metastable  state 
over  a  considerable  temperature  range  above  its  boiling  point,  and  that,  as  a 
result,  the  loss  of  mass  or  energy  by  evaporation  is  negligible  in  picosecond 
laser  quenching. 

In  the  third  part  of  the  paper  the  thermodynamic  restrictions  on  the  type 
of  transformations  that  can  occur  on  the  very  short  timescale  of  the  laser 
quenching  process  are  discussed.  The  selection  of  the  actual  transformation 
from  the  remaining  allowed  ones  depends  on  the  kinetics  of  the  competing 
crystallization  processes.  Depending  on  the  degree  of  structural  rearrange¬ 
ment  required  to  form  a  crystal  from  the  melt,  the  process  can  be  collision- 
or  diffusion-controlled,  it  will  be  shown  that  only  for  the  simplest  crystal 
structures  or  dilute  solutions  is  growth  collision-controlled  and  fast  enough 
to  be  insuppressible,  even  by  picosecond  laser  quenching. 

In  the  last  part  of  the  paper,  two  applications  are  presented.  First, 

It  is  demonstrated  that  fast  pulsed  laser  Irradiation  is  one  of  the  few  methods 
available  for  overheating  metals  far  above  their  melting  temperature.  As  an 
example  of  the  solidification  products  obtained  by  pulsed  laser  quenching,  a 
set  of  recent  experiments  on  the  Nb-SI  system  Is  reviewed.  It  is  an  illustra¬ 
tion  of  the  kinetic  hierarchy,  discussed  by  Turnbull  in  the  first  paper  of 
these  proceedings,  In  which  the  formation  of  isotropic  systems  and  systems 
with  small  unit  cells  Is  klnetically  favored. 

laser  Processing  Mechanism 

A  number  of  detailed  analyses  have  been  made  of  the  processes  of  energy 
deposition,  heating  and  cooling  in  pulsed  laser  irradiation  of  solid  surfaces 
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Fig.  1  -  Schematic  illustration  of  the  mechanism  of  pulsed  laser  quenching. 

(1-3).  These  three  stages  are  illustrated  in  Figure  1,  with  parameters 
corresponding  to  the  experiments  to  be  described  later. 

During  the  first  stage,  the  laser  energy  is  deposited  in  a  layer  of  thick¬ 
ness  a"1,  where  a  is  the  absorption  coefficient.  For  optical  light  on 
metals,  or  for  UV  light  on  semiconductors,  this  layer  is  a  few  hundred  A 
thick.  Since  the  transfer  of  the  energy  from  the  electrons,  which  interact 
with  the  light,  to  the  lattice  occurs  within  less  than  1  ps  (4),  the  energy 
deposition  process,  which  lasts  for  the  duration  of  the  pulse,  tp,  can  be 
considered  a  purely  thermal  one  (5).  The  thermal  diffusion  length  during 
that  time  is  «,T  *  (2Dfht  )1/f2,  which,  for  metals  is  between  250  and  700  %, 
depending  on  the  thermal  diffusivity  (Dth  *  10  to  10  nr/s).  For  fluences, 
4>,  In  the  range  indicated  in  Figure  1,  a  layer  of  thickness  is  melted  and 
heated  to  several  thousand  degrees. 

Note  that  for  a  pulse  duration  of  less  than  1  ps,  £y  is  less  than  a"\ 
The  thickness  of  the  initial  molten  layer,  and  hence  also  the  thermal 
gradients,  are  then  independent  of  the  pulse  duration.  Femtosecond  pulses, 
therefore,  although  very  useful  for  time-resolved  spectroscopy  of  the  various 
phenomena  associated  with  the  irradiation  process,  do  not  produce  a  higher 
quench  rate.  A  general  discussion  of  the  regimes  associated  with  the  relative 
values  of  lj  and  a"^  can  be  found  in  Bloembergen's  paper  (1). 

Following  the  termination  of  the  irradiation  pulse,  the  heat  from  the 
very  hot  liquid  continues  to  diffuse  Into  the  substrate,  which  results  in 
further  melting,  until  the  temperature  of  the  crystal-liquid  in’-.rface  drops 


to  the  equilibrium  melting  temperature.  For  the  conditions  of  Figure  1,  the 
maximum  melt  depth,  d,  is  about  1000  As  will  be  shown  below,  the  over¬ 
heating  of  the  crystal  during  the  melt-in  process  is  considerable.  The 
thermal  gradients  are  very  high  (10^  K/m),  as  is  the  melt-in  velocity 
(u~10  m/s).  The  duration  of  the  melt-in  stage  of  the  process  Is  therefore 
on  the  order  of  d/u~100  ps. 


The  third  stage  is  that  of  resolidification,  either  by  regrowth  of  the 
underlying  crystal,  by  growth  of  a  new  crystal,  or  by  glass  formation  if  the 
melt  cools  to  its  configurational  freezing  point  before  the  crystallization 
front  passes.  A  number  of  sophisticated  numerical  analyses  have  been  made  of 
the  thermal  quantities  of  this  process,  but  their  order  of  magnitude  can 
easily  be  estimated  from  dimensional  arguments.  The  temperature  scale  is  set 
by  the  melting  temperature,  T  ,  which  is  on  the  order  of  10  K,  The  length 
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scale  is  set  by  the  melt  depth  d~10  m.  The  corresponding  temperature 

10  . 

gradient  is  then  VT  *  T^d'-lO  K/m,  and  the  cooling  rate  T  *  Dth  VT/d  ~ 
1012  K/s.  The  lifetime  of  the  melt  can  be  estimated  as  t  *  T^T  *  10‘9  s, 
which  is  In  agreement  with  the  experimental  values  obtained  from  transient 
reflectivity  measurements  (6-9).  The  velocity  with  which  the  isotherms  move 
toward  the  surface  during  cooling  is  an  important  quantity  to  be  compared  to 
the  crystal  growth  velocity,  and  can  be  estimated  as  u-j-  *  T/VT~100  m/s. 
These  parameters,  together  with  those  for  a  "conventional"  melt  quenching 
process  such  as  melt  spinning,  are  sunmarized  in  Table  I,  using  the  thermal 
data  for  iron. 


Table  I. 

Thermal  Parameters  in  Melt  Quenching 

Laser  Quenching 

Melt  Spinning 

Melt  temperature 

V*! 

103 

103 

Melt  thickness 

d(m) 

10"7 

5xl0’5 

Temperature  gradient 

VT(K/m) 

1010 

2x1 07 

Cooling  rate 

T(K/s) 

1012 

4xl06 

Melt  lifetime 

t(s) 

10'9 

not  applicable 

Isotherm  velocity 

uT(m/s) 

100 

0.2 

Heat-flow  limited 
crystal  growth  velocity 

uh(m/s) 

230 

0.5 

The  Overheated  Liquid  State 

As  discussed  in  the  previous  section,  a  molten  layer  with  a  thickness 
JLj.  and  a  temperature  of  several  thousand  degrees  is  formed  in  during  the 
energy  deposition  stage  of  the  pulsed  laser  processing.  Since  the  boiling 


point,  Tb,  of,  for  example,  liquid  iron  is  3160  K  (at  1  atm),  it  is  Important 
to  consider  whether  the  liquid  can  be  overheated  and  hence  be  preserved  for 
later  quenching,  or  is  lost  by  fast  evaporation  above  Tb- 

1 .  Thermodynamic  Stability  Limit 

That  liquids  can  be  overheated  considerably  above  their  boiling  tempera¬ 
ture  Is  a  well-known  experimental  fact  (10,11)  that  can  be  explained  simply 
by  the  van  der  Waals  equation  of  state: 

(p  +  ^-)  (V-b)  *  RT  (1) 

where  a  and  b  are  constants. 


Fig.  2  -  P-V  diagram  for  iron,  showing  one  Isotherm,  the  liquid 
vapor  coexistence  domain,  aid  the  spinodal  lines.  The 
diagram  is  based  on  a  van  Jer  Waals  model  using  the 
critical  point  parameters  of  ref.  (12). 


Figure  2  shows  an  isotherm,  calculated  from  the  critical  point  in  Fe  (12): 


8a/27Rb  =  6753  K 

(2) 

a/27b2  -  501  MPa 

(3) 

3b  *  43  cm^/mole  . 

(4) 

It  should  be  kept  in  mind  that  these  values  for  pc  and  Tc  are  probably 
lower  limiting  ones.  For  example,  a  critical  temperature  as  high  as  10000  K 
has  been  quoted  for  Fe  (12). 

The  equilibrium  pressure  and  volumes  for  the  liquid-vapor  transition  at 
a  given  temperature  are  found  by  the  Maxwell  equal  area  construction,  as 


shown  by  the  cross-hatched  areas  on  the  p-V  diagram  of  Figure  2.  The  locus 
of  the  equilibrium  points  forms  the  coexistence  curve  (solid  line),  inside 
which  the  equilibrium  state  consists  of  a  mixture  of  vapor  and  liquid. 

As  long  as  the  compressibility  of  the  liquid  is  positive,  however,  ( i . e . , 
the  slope  of  the  isotherm  is  negative)  the  liquid  can  be  overexpanded  into 
the  coexistence  region  up  to  S^,  a  point  first  labeled  by  van  der  Waals  as 
the  spinodal  point.  Similarly,  the  vapor  can  be  overcompressed  down  to  the 
spinodal  point  Sy.  The  locus  of  the  spinodal  points  is  easily  obtained  from 
the  equation  of  state  (1),  by  requiring  dp/dV  *  0,  which  gives: 


and  is  indicated  by  the  dashed  line  on  Figure  2. 

In  the  region  between  the  coexistence  line  and  the  spinodal  line,  the 
liquid  or  vapor  are  metastable.  This  means  that  in  order  to  form  the  vapor 
phase  from  the  liquid  It  is  necessary  to  nucleate  a  bubble,  which  occurs  at  a 
finite  frequency.  Similarly,  in  order  to  form  a  liquid  from  the  metastable 
vapor  requires  nucleation  of  a  droplet.  Inside  the  spinodal  line  the  nuclea- 
tion  barrier  vanishes:  a  fluctuation  of  any  size  is  stable,  and  vaporization 
or  condensation  occurs  rapidly  by  spinal  decomposition  of  the  liquid  or  vapor 
into  two  phases. 

Figure  3  shows  the  metastable  regions  for  the  liquid  and  vapor  of  Fe  on 
a  p-T  diagram.  The  spinodal  lines  and  the  upper  part  of  the  equilibrium  line 
(solid  line)  have  been  calculated  from  the  van  der  Waals  model  and  the 
critical  parameters  listed  above.  The  dashed  part  of  the  equilibrium  line 
has  been  adjusted  to  go  through  the  boiling  point  at  1  atm  (0.1  MPa). 


In  laser  heating  of  a  molten  surface  at  atmospheric  pressure,  rapid 
vaporization  by  spinodal  decomposition  throughout  the  liquid  can  only  occur 
if  the  temperature  exceeds  that  of  the  liquid  spinodal  line  at  1  atm,  T$ 

This  can  be  calculated  precisely  by  putting  p*0  in  equation  (5),  and  using 
the  critical  point  relations  (2)-(4): 


(6) 


For  iron,  using  the  lower  limiting  critical  parameters  (2 )- (4 ) ,  this 
temperature  is  estimated  to  be  at  least  5700  K,  as  can  be  checked  on  Figure  3. 
It  could  actually  be  as  high  as  8400  K,  if  Tc  is  10000  K. 

If  the  temperature  of  the  liquid  exceeds  T$  ^ ,  the  decomposition  pro¬ 
ceeds  very  rapidly,  but  an  analysis  of  Its  time  dependence  requires  a  quanti¬ 
tative  kinetic  theory.  The  theory  of  spinodal  decomposition  is  well 


Temperature  (K5 

Fig.  3  -  P-T  diagram  showing  the  liquid-vapor  equilibrium  line, 

the  spinodal  lines  and  the  metastable  regions  for  liquid 
and  vapor.  The  diagram  is  based  on  a  van  der  Waals 
model  using  the  critical  point  parameters  of  ref.  (12). 

developed  (13-15),  and  explains  the  minimum  length  scale  for  decomposition  by 

taking  into  account  the  contribution  of  very  steep  gradients  (in  density  or 

concentration)  to  the  free  energy  of  an  inhomogeneous  system.  Although  an 

estimate  of  the  gradient  energy  coefficient  for  liquid  metals  can  in  principle 

be  obtained  from  the  same  mean  field  interactions  that  lead  to  the  van  der 

Waals  model,  we  can  simply  use  the  average  interatomic  distance  in  the  equi- 

0 

librium  vapor  phase  at  T_  0  (~6  A)  as  a  lower  limit  for  the  length  scale  of 
the  density  fluctuations.  An  approximate  argument  similar  to  the  one  used 
below  to  estimate  the  time  required  for  assembling  a  critical  vapor  nucleus 
(section  3)  can  then  be  invoked  here  for  the  formation  of  a  fluctuation.  It 
is  therefore  possible  that,  because  of  the  very  short  time  scale  of  picosecond 
pulsed  laser  irradiation,  the  liquid  can  even  sustain  some  overheating  above 


2.  Steady  State  Nucleation 

If  the  temperature  of  the  molten  layer  does  not  reach  T  formation 

5  f  Ar 

of  the  vapor  Inside  the  liquid  requires  nucleation  of  a  bubble.  Classical 
nucleation  theory  (1S-18)  give  the  following  expression  for  the  eteady  state 
nucleation  frequency  (number  of  nuclei  per  unit  volume  and  time): 


where  N  Is  the  number  of  atoms  per  unit  volume,  k  is  Boltzmann's  constant, 
h  is  Planck's  constant,  p  is  the  equilibrium  pressure  at  the  temperature  T 
(p  is  much  larger  than  the  ambient  pressure  here),  and  o  is  the  liquid- 
vapor  surface  tension. 

The  minimum  nucleation  frequency  required  to  cause  vaporization  through¬ 
out  the  liquid  layer  in  laser  quenching  can  be  estimated  as: 


min 


(8) 


where  a  is  the  maximum  distance  traveled  in  a  time  tp  by  the  vapor-liquid 

interface  of  a  nucleated  bubble;  a  =  1  ym  is  taken  as  an  absolute  upper 

limit.  This  gives  a  value  for  I  .  of  2*1030  m‘3s'1.  The  maximum 

temperature  to  be  considered  here  is  Tc  c,  so  that  an  estimated  T =  5500  K 

is  used  in  the  analysis  below;  it  should  be  kept  in  mind  that  the  results  are 

very  insensitive  to  the  choice  of  the  values  for  T  and  I  .  .  The  equili- 

mi  n 

brium  vapor  pressure  at  this  temperature  can  be  read  from  Figure  3: 
p  «=  200  MPa. 

Equation  (1)  can  now  be  used  to  calculate  a  maximum  value  for  the  liquid 
vapor  surface  tension,  o,  for  the  nucleation  frequency  to  exceed  Imin-  This 
gives  a  value  o  *  0.17  N  m’\  The  surface  tension  of  liquid  iron  at  the 

TTlaX  ^ 

melting  point  (1812  K)  is  1.788  N  m"  (19).  Since  the  surface  tension 
decreases  roughly  linearly  with  increasing  temperature  and  vanishes  at  the 
spinodal  T_  „  (14),  the  surface  tension  drops  below  om.v  only  if  the 
temperature  comes  within  a  few  hundred  degrees  from  T  .. 


3.  Transient  Nucleation 

Equation  (7)  describes  the  isothermal  steady  state  nucleation  frequency, 
l.e.,  the  rate,  per  unit  volume,  at  which  bubbles  form  if  a  steady  state 
population  of  vapor  embryos  has  been  established.  Establishing  such  a  popula 
tion  after  a  sudden  change  in  temperature,  however,  requires  a  certain  amount 
of  time,  during  which  the  nucleation  rate  increases  to  its  steady  state 


value  (20).  Usually,  this  transient  effect  is  only  important  in  transforma¬ 
tions  In  condensed  systems,  where  the  atomic  mobility  is  relatively  low.  In 
vapor-liquid  transitions,  the  transient  time  is  much  shorter  and  can  usually 
be  ignored.  On  the  very  short  time  scale  of  pulsed  laser  irradiation,  how¬ 
ever,  transient  effects  may  be  important,  and  may  be  a  stabilizing  factor 
against  copious  nucleation  in  the  overheated  liquid. 

Although  an  exact  analysis  of  the  transient  problem  is  very  time- 
consuming,  a  lower  limit  of  the  transient  time,  t,  can  be  estimated  by  con¬ 
sidering  the  time  it  takes  to  form  a  critical  vapor  nucleus  by  evaporation  of 
the  liquid. 

For  large  equilibrium  vapor  pressures  p,  the  critical  nucleus  radius  is 
given  in  the  classical  theory  by: 


which,  for  the  nucleation  example  above  is  17  %. 


The  rate  of  evaporation  can  be  estimated  from  the  kinetic  theory  of 
gases.  The  maximum  number  of  atoms  evaporating  from  a  unit  surface  per  second 
is,  in  equilibrium,  equal  to  the  number  impinging  on  that  surface  from  the 
vapor,  which  is: 


r  *  — -P — 

/SirinkT 


(10) 


where  m  is  the  atomic  mass.  For  the  above  conditions,  this  is  about 
30  -2  -1 

10  m  s  .  If  the  sticking  coefficient,  q,  of  the  vapor  atoms  is  less  than 
unity,  T  is  reduced  to  (17): 

r'  ■  r  2?q  •  (11) 

The  time  to  form  a  spherical  nucleus  of  radius  r*  can  then  be  estimated 
by  considering  that  at  any  time  during  the  growth  of  the  vapor  bubble  the 
total  volume  of  the  atoms  evaporating  from  the  surface  of  the  bubble  (of 
Instantaneous  area  A)  is  equal  to  the  instantaneous  increase  in  the  bubble 
volume,  or: 

Afirdt  *  A  dr  (12) 

where  fl  is  the  atomic  volume  of  the  atoms  in  the  vapor  phase.  Figure  2 
shows  that  for  the  conditions  considered  here  the  molar  volume  of  the  vapor 
Is  around  150  cm^,  which  corresponds  to  8*  250  ft.  Equation  (12)  then 
integrates  simply  to  give  the  minimum  time  to  form  the  critical  nucleus: 

rl 

8T 


T  * 
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The  holes,  therefore,  are  not  a  result  of  evaporation  per  ee ,  but  of 
mechanical  displacement  of  the  liquid  layer  by  the  large  recoil  pressure  of 
the  evaporation  ("splashing").  The  morphology  of  the  displaced  material 
around  the  hole  is  often  direct  evidence  of  this  splashing  phenomenon,  as  is 
the  finding  of  substrate  atoms  near  the  hole  on  the  top  side  of  supported 
thin  films  (21). 

Direct  evidence  of  the  recoil  pressure  can  be  obtained  from  the  following 
experiment  (22).  A  thin  ribbon,  clamped  at  one  end,  is  irradiated  by  a  laser 
pulse  at  the  free  end,  which  results  In  visible  vibration.  The  impulse 
(force  x  time)  of  the  recoil  pressure  calculated  from  the  vibration  amplitude 
is  in  general  agreement  with  the  order  of  magnitude  of  the  temperatures  and 
times  discussed  above. 


Thermodynamics  and  Kinetics 

The  short  time  scale  of  the  picosecond  melt  quenching  process  puts  re¬ 
strictions  on  the  transformations  that  can  occur.  For  the  entire  liquid  layer, 
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of  thickness,  d*1000  A,  to  crystallize  during  the  cooling  time,  t,  of  1  ns, 
the  velocity  must  be  at  least,  u  *  d/x  ~  100  m/s.  The  time,  t^,  to 
crystallize  a  monolayer  at  this  speed  is  about  3  ps,  and  the  distance  an  atom 
can  diffuse  in  that  time,  (D^t-j)  '  ,  is  less  than  an  interatomic  distance. 

(0£  is  the  diffusivity  in  the  liquid  phase).  This  means  that  under  these 
short  pulsed  laserquenching  conditions  no  long-range  diffusion  can  occur,  and 
that  the  only  possible  transformations  are  partitionless  ones,  either  growth 
of  a  crystal  of  the  same  composition  of  the  liquid,  or  glass  formation. 

The  thermodynamic  restrictions  imposed  by  the  partitionless  conditions 
are  Illustrated  by  the  phase  diagram  of  Figure  4,  which  contains  two  primary 
solid  solutions  (a, 6)  and  two  intermetallic  compounds  (0,y).  The  free 
energy  diagram  In  Figure  4  shows  how  the  TQ-lines  for  the  transformations  of 
these  four  phases  from  the  liquid  are  constructed.  At  a  sufficiently  low 
temperature,  Tg,  the  atomic  transport  rate  in  the  liquid  state  becomes 
negligibly  small,  so  that  no  more  transformations  can  occur  on  a  reasonable 
time  scale.  The  liquid  Is  then  configurationally  frozen  and  becomes  a  glass. 
For  a  composition  In  the  range  x'  -  x"  in  Figure  4,  no  crystallization  is 
possible  under  picosecond  laser  quenching  conditions,  and  therefore  only  glass 
formation  is  thermodynamically  possible.  This  criterion  for  glass  formation 
Is  only  a  sufficient  one,  and  not  a  necessary  one,  as  has  been  claimed  (23,24). 
In  fact,  more  often  than  not,  picosecond  quenching  experiments  yield  glasses 
over  a  wide  range  of  compositions  below  the  T0-lines  for  both  intermetallic 
compounds  and  primary  solutions.  To  explain  this,  it  is  necessary  to  consider 
the  kinetics  of  partitionless  crystallization  in  some  detail. 


Fig.  4  -  Schematic  phase  diagram  (top),  and  corresponding  free 

energy  diagram  at  temperature  T]  (bottom),  illustrating 
the  construction  of  the  T^-lines  for  the  primary  solid 
solutions  and  the  intermetall ic  compounds.  In  the  compo¬ 
sition  range  x‘ -  x"  (dashed),  only  glass  formation  is 
possible  under  conditions  of  partitionless  solidification. 
The  numbers  correspond  to  the  same  compositions  in  the 
top  and  bottom  diagrams. 


The  basic  theory  for  crystal  growth  far  from  equilibrium  has  been  formu¬ 
lated  by  Turnbull,  and  can  be  found  in  a  number  of  publications  (25,26), 
including  the  first  paper  in  the  present  volume.  The  crystal  growth  velocity 
can  be  written  as: 


where  k  is  the  atom  jump  frequency  across  the  crystal -melt  interface,  f  is 
the  fraction  of  interface  sitts  that  can  incorporate  a  new  atom,  X  is  the 
Interatomic  distance,  &GC  is  the  difference  in  molar  free  energy  between 
crystal  and  melt,  and  T^  is  the  interface  temperature. 

If  the  undercooling  of  the  interface,  =  Tm  -  T^ ,  is  small,  the 
driving  free  energy  can  be  approximated  by: 
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AGc  *  AScATi 

where  AS£  Is  the  molar  entropy  of  crystallization.  The  exponential  in 
equation  (15)  can  be  linearized  to: 


u  *  fkX 


ASC  AT, 


(17) 


where,  for  metals,  f~l  and  AS£~R. 

The  crystal  growth  rate  also  depends  on  the  rate  at  which  the  heat  of 
crystallization,  Ah  per  unit  volume,  can  be  removed.  Since  this  must 

C  »  V 

occur  by  heat  flow  down  the  temperature  gradient  at  the  interface,  V^T,  the 
heat  flux,  J,  must  obey  the  relations: 

J  *  uAh  w  «  kV.T  (18) 

t  »  V  1 


where  k  is  the  thermal  conductivity.  Conversely,  equation  (18)  can  be 
turned  around  to  define  a  heat  flow-limited  velocity: 


(19) 


Table  I  gives  values  for  uh  for  the  laser  quenching  conditions  as  well  as 
for  melt  spinning. 

Oomoining  equations  (17)  ^d  (19)  then  gives  for  metals: 


VTi 


(20) 


In  conventional  solidification,  the  interface  kinetics  are  much  faster  than 
the  rate  of  heat  removal,  so  that  u^  «  kX,  and  AT^  is  small.  This  type 
of  growth  is  called  heat  flow-limited.  In  laser  quenching,  depending  on  the 
nature  of  the  interface  rearrangements  (see  below),  the  rate  of  heat  removal 
can  be  much  faster  than  the  interface  kinetics,  so  that  uh  »  kX,  and  the 
undercooling,  AT^,  is  large,  so  that  metastable  phases  can  be  formed.  This 
corresponds  to  interface-limited  growth. 

The  maximum  growth  velocity  of  a  crystal,  u  s  kX,  is  determined  by 

the  nature  of  the  atomic  rearrangements  that  control  the  growth  process.  In 

pure  metals,  dilute  alloys  or  compounds  with  simple  crystal  structures,  k 

can  be  taken  as  the  frequency  of  thermal  vibration.  (~  Debye  frequency),  so 

that  u„„w  ~  u  ,  the  speed  of  sound  in  the  liquid,  which  is  on  the  order  of 
max  s  r 

several  thousand  m/s.  This  Is  called  the  collision- controlled  regime, 
because  the  growth  occurs  by  simple  impingement  of  the  liquid  ato  3  onto  the 
crystal . 


In  more  concentrated  alloys  or  in  compounds  with  a  more  complicated 
crystal  structure,  where  a  change  of  the  nearest  neighbor  configuration,  i.e. 
diffusion,  is  necessary,  the  jump  frequency  can  be  taken  as  D^/X  ,  and 
Umax  d  *  Di/*'  which  on  th€  order  of  10  m/s.  This  is  the  diffusion - 
controlled  regime.  Table  I  shows  that  for  picosecond  laser  quenching 

umax  D  <<  UT’  uh*  so  that  this  type  9rowth  can  easily  be  suppressed. 

Equation  (20)  shows  that  under  these  conditions  T^  falls  far  below  Tm,  and 
glass  formation  becomes  possible. 

Although  in  partitionless  growth  of  crystals  there  is  no  long-range 
diffusional  transport,  the  growth  can  still  be  diffusion-controlled  if  a 
drastic  change  of  the  nearest  neighbor  environment  is  necessary  to  transform 
the  liquid  structure  into  the  crystalline  one.  This  is  the  case,  for  example, 
in  metal-metalloid  alloys  such  as  Fe-B,  where  the  chemical  interaction  leads 
to  a  very  pronounced  short-range  order  in  the  liquid  (9  Fe  atoms  around  each 
B  In  a  capped  trigonal  prism  arrangement  r27)),  which  is  very  different  from 
that  in  the  crystal  (B  interstitial  in  bcc  Fe  (28)).  As  a  result,  it  is 
possible  to  produce  Fe-B  glasses  containing  as  little  as  5  at.%  B  (28),  which 
is  far  below  the  TQ-11ne.  Similarly,  if  the  unit  cell  of  an  Intermetal  1 ic 
compound  is  large,  or  if  long-range  chemical  order  must  be  established,  crystal 
growth  is  diffusion-limited,  and  the  growth  of  these  compounds  is  easily 
suppressed  in  picosecond  pulsed  laser  quenching.  Numerous  experiments  on  a 
variety  of  binary  metallic  systems  (Ni-Nb,  Ni-Mo,  Co-Mo,  Nb-Co  (8,  30-34)) 
have  demonstrated  that  glasses  can  be  formed  at  the  intermetallic  compositions. 
On  the  other  hand,  in  pure  metals  with  simple  crystal  structures,  and  by 
extension,  disordered  alloys  with  the  same  simple  structures  (e . g. ,  disordered 
fee  in  Cu-Co  or  Ag-Cu),  growth  Is  collision-limited  and  cannot  be  suppressed 
even  by  picosecond  pulsed  laser  quenching.  This  Is  confirmed  by  direct 
measurements  of  the  growth  velocity  of  pure  metals  and  disordered  alloys  (Au, 
Cu,  Cu-Au)  by  transient  reflectivity  (8,9)  which  yield  large  values  (>  100  m/s) 
that  can  only  be  accounted  for  by  a  collision-controlled  mechanism. 

Appl ications 


1 .  Overheating  of  Crystals 

It  is  well  known  that  the  interior  of  large  crystals  can  be  heated  well 
above  the  melting  temperature  if  the  kinetics  of  melting  are  much  slower  than 
the  rate  of  thermal  diffusion,  as  Is  the  case  for  oxides  (35).  In  metals  the 
kinetics  of  melting  are  much  faster,  so  that  overheating  Is  much  more  difficult 
to  achieve.  One  method.  In  which  conventional  heating  is  used,  consists  of 
surrounding  a  small  single  crystal  with  a  higher  melting  coherent  layer  to 


prevent  surface  nucleation  of  the  interior  melt.  Overheating  of  silver 
crystals,  coated  by  a  gold  layer,  by  7.5  K  has  been  reported  (36).  ‘Conven¬ 
tional"  rapid  heating  methods,  for  example  by  sending  a  high  current  pulse 
through  a  thin  metallic  wire,  also  give  overheating  of  no  more  than  a  few 
degrees  (37). 

Rapid  heating  with  pulsed  lasers  provides  a  unique  method  for  overheating 
crystals  considerably  auove  their  melting  temperature,  as  can  be  seen  from  the 
following  argument.  During  the  Initial  stages  of  the  melt-in  phase  the 
temperature  gradients  are  an  order  of  magnitude  greater  than  during  regrowth. 
This  means  that  the  heat-flow  limited  crystal  growth  velocity,  given  by 
equation  (19),  Is  also  an  order  of  magnitude  greater  than  the  value  listed  in 
Table  I,  or  about  2000  m/s.  For  pure  metals,  where  u„v  «  ue  ~  3000  m/s, 
equation  (20)  gives  for  the  overheating:  ■  0.66  T^,  which,  for  a  melting 

point  of  1000  K,  corresponds  to  an  overheating  of  2000  K. 

The  obvious  drawback  of  the  method  is  the  short  duration  of  the  over¬ 
heating  (<  100  ps),  which  requires  fast  diagnostic  techniques  for  studying 
the  crystal  or  possible  homogeneous  nucleation  of  the  melt  in  this  state. 

2.  Quenchino  of  Nb-Si  Alloys 

The  Nb-Si  system  is  of  interest  because  of  the  prediction  that  the  Nb^Si 
compound  with  the  A-15  structure  should  have  a  very  high  superconducting 
transition  temperature.  Since  at  25  at.%  Si  there  are  no  stable  compounds  in 
this  system  (see  the  phase  diagram  of  Figure  5),  a  variety  of  methods  have 
been  used  to  produce  metastable  phases.  As  Illustrated  in  Figure  5,  conven¬ 
tional  melt  quenching  methods  produce  an  amorphous  phase  in  the  vicinity  of 
the  eutectics  (38).  Non-stoichiometric  A-15  phases  have  also  been  found 
(39,40).  At  the  25  at.%  Si  composition,  however,  conventional  melt  quenching 
yields  a  compound  with  the  tetragonal  TigP  structure  (41). 

We  have  recently  completed  a  study  of  this  system  using  pulsed  laser 
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quenching  of  thin  alloy  films  (~  1000  A)  supported  on  various  substrates 
(42-44).  Although  we  did  not  succeed  In  obtaining  an  A-15  phase  (the  quench 
rates  were  too  high  in  all  attempts  for  this),  it  Is  still  useful  to  review 
the  various  phases  that  were  obtained  at  different  cooling  rates.  Figure  5 
Illustrates  the  results. 

1 2 

At  the  highest  quench  rates  (>  10  K/s),  obtained  in  the  thinnest  sec¬ 

tions  of  the  spots  irradiated  with  a  30  ps  pulse  on  a  copper  substrate  the 
amorphous  phase  was  formed  over  a  wide  composition  range.  Such  broadening  of 
the  range  of  amorphous  phase  formation  is  consistent  with  our  earlier  observa¬ 
tions  on  other  metal -metalloid  and  metal -metal  systems. 
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Fig.  5  -  Phase  diagram  of  the  Nb-Si  system  (top),  and  schematic 
illustration  of  the  phases  obtained  by  different  methods 
of  liquid  quenching.  The  melt  spinning  results  are  from 
refs.  (38-41);  the  laser  quenching  results  are  from 
refs.  (42-44).  The  pulse  length  and  the  substrate  for 
the  laser  experiments  are  indicated. 

1 2 

At  slightly  lower  quench  rates  (10  K/s),  obtained  in  the  remaining 

area  of  those  spots,  a  new  phase  was  found:  a  disordered  f.c.c.  solid  solu¬ 
tion  (43).  Although  the  formation  of  disordered  f.c.c.  solid  solutions  in 
metal-metal  alloys  cannot  be  suppressed  by  picosecond  laser  quenching,  the 
growth  rate  of  the  f.c.c.  phase  in  this  system,  although  still  the  highest 
one  observed,  is  probably  slightly  lower  due  to  the  short-range  order  around 
the  Si. 
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By  lowering  the  quench  rate  to  10  -  10  K/s  by  supporting  the  films  on 

an  aluminum  substrate,  a  supersaturated  b.c.c.  solid  solution  was  found  at 
low  SI  contents,  and  an  ordered  f.c.c.  phase  (L^,  Cu^Au-type)  around  25  at.% 
Si. 

9  10 

A  cooling  rate  of  10  -  10  K/s  was  obtained,  either  by  supporting  the 
ps-irradlated  films  on  a  salt  substrate,  or  by  Irradiating  the  films  on  a 
copper  substrate  with  a  30  ns  pulse.  Under  these  conditions  two  entirely  new 
phases  were  found  around  25  at.%  Si:  a  hexagonal  phase  with  a  250  unit 


cell  and  an  as  yet  unidentified  phase  (44). 

It  is  interesting  to  note  the  degree  of  isotropy  and  the  unit  cell  size 
of  the  phases  formed  at  25  at. Si  SI  in  order  of  decreasing  quench  rate: 
amorphous  (isotropic,  no  unit  cell),  disordered  f.c.c.  (1  atom  unit  cell, 
cubic  point  group),  ordered  f.c.c.  (4  atom  unit  cell,  cubic  point  group), 
hexagonal  (~  16  atom  unit  cell,  uniaxial  point  group),  Tl^P-type  (32  atom 
unit  cell,  uniaxial  point  group).  This  is  a  good  illustration  of  the  prin¬ 
ciple,  stated  in  Turnbull's  paper  in  this  volume,  that  phases  are  increasingly 
klnetically  favored  (i.e.,  are  more  difficult  to  suppress  by  quenching)  with 
Increasing  Isotropy  and  decreasing  unit  cell  size. 
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Crystallization  rates  in  pure  metals  have  been  measured  by  monitoring  changes  in  reflectivity 
or.  a  -iccsecond  time  scale  during  pulsed  hser  quenching  Multilayered  films  were  used  to 
establish  that  the  transient  decreases  in  reflectivity  correspond  to  melting.  Depth  profiling  with 
Auger  electron  spectroscopy  on  laser-irradiated  areas  showed  complete  mixing  of  the  layers, 
which  could  only  have  occurred  by  melting.  Crystallization  velocities  as  high  as  100  m/s  were 
measured.  These  growth  speeds  are  too  large  to  be  the  result  of  a  crystallization  process 
governed  by  liquid  state  diffusion.  The  results  of  computer  simulations  based  on  a  collision- 
limited  growth  model  for  crystallization  were  consistent  with  the  measurements. 


I.  INTRODUCTION 

Short  pulse  lasers  (tens  of  picoseconds)  deposit  heat  in  a 
very  thin  region  near  the  surface  of  a  sample.1,2  If  sufficient 
energy  is  deposited,  a  melt  nucleates  at  the  surface  and  pro¬ 
pagates  inward.  During  this  period  the  underlying  crystal¬ 
line  material  can  be  superheated  substantially  above  the 
melting  temperature.  After  the  end  of  the  laser  pulse  the 
molten  layer  begins  to  cool  as  heat  flows  into  the  interior  of 
the  sample.  When  the  temperature  of  the  melt/crystal  inter¬ 
face  drops  below  the  melting  temperature,  the  interface  re¬ 
verses  direction  and  returns  to  the  surface.  The  large  thermal 
gradients  present  in  short  pulsed  laser  quenching  experi¬ 
ments  can  result  in  substantial  undercooling  of  the  interface 
below  the  melting  temperature.  A  complete  theory  for  the 
kinetics  of  crystallization  that  adequately  predicts  the  veloc¬ 
ity  of  this  interface  as  a  function  of  its  temperature  is  still 
being  developed.3-5  There  has  been  extensive  work  devoted 
to  measurements  of  interface  velocities  in  silicon  by  monitor¬ 
ing  transient  optical  and  electrical  properties  during  and 
after  the  laser  pulse.6  These  experiments  were  facilitated  by 
the  dissimilarity  of  the  crystal  and  its  melt.  Upon  melting, 
silicon  changes  from  a  covalently  bonded  semiconducting 
solid  with  four  nearest  neighbors  to  a  metallic  liquid  with 
high  coordination  number.  The  normal  incidence  reflectiv¬ 
ity  differs  by  a  factor  of  2  and  the  dc  conductivity  by  a  factor 
of  30.  On  the  other  hand,  the  theoretical  analysis  of  crystalli¬ 
zation  in  silicon  is  hampered  by  this  same  dissimilarity. 
Further,  since  the  structural  change  upon  crystallization  is 
simpler  in  metals  than  in  silicon,  it  might  be  expected  to  be 
faster.  In  order  to  achieve  high  velocities,  it  is  necessary  to 
produce  large  thermal  gradients  such  as  those  that  occur 
following  picosecond  pulsed  laser  quenching.  This  work, 
and  previous  measurements  by  this  group,7,8  constitute  the 
first  measurements  of  the  kinetics  of  crystallization  in  metals 
during  pulsed  laser  quenching,  and  the  only  ones  to  date  in 
the  picosecond  regime. 

•’Currently  at  Physics  Dept  ,  Tne University  at  Albany,  State  University  of 
New  York,  Albany,  NY  12222. 

*”  Currently  at  Dipartimentodi  Elettronica,  Universitadi  Pavia.  Via  Abbia- 
tegrasso  209,  Pavia,  Italy 


II.  EXPERIMENTAL  TECHNIQUE 
A.  Optical 

The  experiment  was  conducted  with  a  standard  pump/ 
probe  optical  setup,  shown  in  Fig.  1.  The  30-ps  Nd:YAG 
pulse  at  1.06/zm  was  doubled  to  produce  a  20-ps  green  (0.53 
//m)  pump  pulse,  used  to  melt  the  sample.  The  spot  sizes 
were  a  few  hundred  microns  in  diameter.  The  reflectivity  of 
the  sample  was  probed  in  the  center  of  the  irradiated  area 
with  a  small  fraction  of  the  remainder  of  the  1  06-/jm  pulse  at 
a  variable  delay  relative  to  the  pump  pulse.  Data  from  this 
system  were  collected  with  an  IBM  PC  XT.9  By  closing  a 
shutter  to  block  the  pump  beam,  the  probe  beam  could  be 
fired  alone  to  determine  the  initial  and  final  reflectivity  of 
each  spot.  A  fresh  portion  of  the  sample  was  irradiated  at 
each  measurement.  To  reduce  the  effect  of  spot-to-spot  var¬ 
iations  in  reflectivity,  the  measurements  were  normalized  to 
the  initial  value  for  each  spot.  Figure  2  displays  a  typical  plot 
of  the  change  in  reflectivity  with  a  delay  of  2  ns  (crosses)  and 
of  several  seconds  (circles)  between  the  pump  beam  and 
probe  pulses.  For  low  fluences  (region  I),  there  was  no 
change  in  reflectivity.  For  moderate  fluences  (region  II),  a 
small  change  in  reflectivity  was  observed  at  2  ns,  with  nearly 


VIG  1  Schematic  diagram  of  the  lime-resolved  reflectivity  apparatus  Dl- 
D5  detectors;  DC:  doubling  crystal 


129 


J  Aool  Phys  65(1).  1  January  1989 


0021-8979/89/0101 29-08S02  40 


<c  1 988  American  Institute  of  Physics 


129 


0.2 


AR 

Ro 


-0.5 


-0.8 


1 

1 

*  2000p$  deloy 
of  mol 

i  •  v  ♦ 

6°°CP  O 

< 

1 

'  1 

*♦** 

-ill 

1 

1 

I  ♦ 

_ ii, 

■  -  1 _ _ _ L  _ 

0.1  0.2 
Fluence  (J/cmJ) 


0.3 


FIG.  2  Relative  change  in  reflectivity  of  a  copper/gold  artificial  multilayer 
with  10-nm  repeat  length,  irradiated  w  ith  a  20- ps  pulse.  Crosses:  2-ns  delay; 
Circles:  several  seconds  delay. 


complete  recovery.  Finally,  for  large  fluences  (region  III), 
large  and  permanent  drops  in  reflectivity  were  observed. 

To  analyze  this  data  it  was  necessary  to  determine 
whether  melting  leads  to  a  decrease  or  increase  in  sample 
reflectivity.  Simple  free-electron  Drude  calculations  predict 
that  the  reflectivity  is  unity  below  the  plasma  frequency,  and 
zero  above.  In  a  real  (finite  temperature)  system,  electron- 
phonon  collisions  reduce  the  reflectivity  for  frequencies  be¬ 
low  the  plasma  frequency.  Since  the  plasma  frequency  for 
metals  is  in  the  near  UV  (or  in  the  blue-green  region  for  the 
yellow  metals  like  copper  and  gold),  the  probe  frequency 
( 1.06/zm)  is  well  below  the  plasma  frequency  for  these  ex¬ 
periments.  Increases  in  temperature  further  reduce  the  re¬ 
flectivity.  On  the  other  hand,  melting,  with  its  associated 
disorder,  increases  the  electron-phonon  interaction,  and 
hence  is  also  expected  to  lower  the  reflectivity.  Thus,  the 
changes  in  reflectivity  seen  in  the  middle  region  in  Fig.  2 
might  a  priori  be  associated  with  either  heating  or  melting. 
Since  small  changes  in  film  preparation  can  result  in  large 
differences  in  reflectivity,  the  literature  contains  widely 
varying  values  for  reflectivity  changes  upon  melting,  which 
are  therfore  not  useful  for  direct  comparison  with  the  pres¬ 
ent  experimental  data.10 


B.  Sample  preparation  and  characterization 

The  experiment  was  performed  on  bulk  specimens  of 
copper,  aluminum,  and  silver,  and  on  thin  films  of  copper, 
gold,  silver,  and  aluminum  on  a  variety  of  substrates  (in 
order  of  increasing  conductivity:  glass  slides,  polished  sili¬ 
con  wafers,  and  sapphire  crystals).  All  samples  exhibited 
similar  behavior.  The  bulk  samples,  while  providing  faster 
cooling  rates  because  of  their  higher  thermal  conductivities, 
exhibited  larger  spot-to-spot  variations  in  reflectivity  and 
localized  melting  on  polishing  scratches  or  etch  pits  left  from 
electropolishing.  Films  produced  by  vapor  deposition  in  a 
chamber  with  no  facility  for  in  situ  cleaning  tended  to  adhere 
poorly  to  the  substrate.  Better  results  were  obtained  by  rf 
sputtering  in  a  SEM  coatcr  after  in  situ  plasma  etching.  Fig¬ 
ure  3  displays  a  typical  100-nm-thick  gold  film,  deposited  by 
the  latter  method  on  a  silicon  substrate,  after  laser  irradia¬ 
tion.  The  ring  pattern  resulted  from  the  circularly  symmet¬ 
ric  Gaussian  intensity  profile  of  the  beam 

The  best  samples  (with  high  reflectivity,  good  adhesion, 
and  low  spot-to-spot  variation)  were  obtained  by  ion  beam 


FIG.  3  Optical  micrograph  of  a  0. 1-^m-thick  gold  film  on  a  silicon  sub¬ 
strate  irradiated  with  a  20-ps  pulse,  (a)  Ring  pattern  resulting  from  the 
Gaussian  intensity  distribution  of  the  beam:  (b)  detail  of  the  laser  damage 
at  the  center  of  the  spot. 


sputtering  in  a  previously  described  system."  This  system 
not  only  allowed  deposition  of  the  largest  variety  of  materi¬ 
als,  but  also  the  production  of  artificial  multilayered  films. 
These  films  were  used  to  determine  if  the  characteristic  re¬ 
flectivity  decreases  could  be  associated  with  melting.  Alter¬ 
nate  layers  of  materials  that  are  completely  miscible  in  the 
liquid  and  solid  states  (copper/gold  and  copper/nickel) 
were  produced  with  repeat  lengths  of  2.5  to  20  nm,  which 
were  measured  by  x-ray  diffraction. 13  Solid-state  diffusion 
rates  are  too  slow  to  allow  substantial  mixing  during  the 
short  time  in  which  the  films  were  hot  after  a  laser  pulse.  The 
diffusivity  in  these  metals  just  below  their  melting  tempera¬ 
ture  does  not  exceed  ~10_l<cm:/s.  If  the  sample  remains  at 
this  temperature  for  1  ns,  the  diffusion  distance  (£>f)l/:  is 
only  0.03  nm.  Thus,  if  substantial  mixing  of  the  film  is  ob¬ 
served  after  laser  irradiation,  melting  must  have  occurred.  A 
comparison  of  the  Auger  electron  spectroscopy  (AES)  sput¬ 
ter  profiles  of  a  copper/nickel  multilayer  of  10  nm  repeat 
length  near  and  inside  a  laser  spot  ( Fig.  4 )  clearly  shows  the 
mixing.  In  order  to  make  valid  comparisons  between  the 
results  from  multilayered  and  pure  films,  it  must  be  estab¬ 
lished  that  they  have  similar  optical  properties.  This  was 
done  by  calculating  the  reflectance  of  the  multiple  layer 
structure  shown  in  Fig.  5.  first  for  a  pure  gold  film,  and  then 
fo:  a  metal  film  consisting  of  alternate  layers  of  gold  and 
copper.  The  coefficients  of  reflectance  and  the  forward  and 
back  transmission  were  computed"  from  the  values 
n  =  0.194,  k  =  5.57  for  solid  gold"  and  n  =  1.38,  k  =  7.2 
for  liquid  gold,"  and  were  used  recursively  to  calculate  the 
overall  reflectance  of  the  system  as  a  function  of  the  depth  of 
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FIG.  4.  Auger  electron  spectroscopy  depth  profile  of  the  copper/gold  mul¬ 
tilayer  of  Fig.  2  The  solid  line  is  for  an  area  adjacent  to  a  laser  spot.  The 
dotted  line  is  for  an  area  irradiated  with  the  fluence  marked  in  Fig.  2. 


the  molten  layer.  The  results  for  pure  gold  are  shown  in  Fig. 
6(a).  The  reflectivity  falls  offby  a  few  percent,  and  levels  off 
at  a  liquid  depth  equal  to  the  optical  absorption  length.  Fig¬ 
ure  6(b)  displays  the  results  for  goldAopper  multilayer  of 
20-nm  repeat  length,  using  n  =  0.92,  k  =  7.7  for  solid  cop¬ 
per  and  n  =  0.43,  k  —  5.6  for  liquid  copper  (from  the  same 
references  as  the  respective  gold  data). 

Figure  2  shows  the  reflectivity  change  of  a  copper/gold 
multilayer  after  a  laser  pulse  as  a  function  of  fluence.  This 
reflectivity  behavior  was  also  typical  of  pure  films  and  bulk 
specimens.  An  AES  profile  of  a  spot  irradiated  with  the 
fluence  marked  in  Fig.  2  (in  the  center  of  region  II)  is  shown 
in  Fig.  4.  The  AES  profile  clearly  shows  that  region  II  can  be 
associated  with  melting.  Bright-  and  dark-field  optical  mi¬ 
crographs  of  the  spot  profiled  in  Fig.  4  are  shown  in  Figs. 
7(a)and7(b).  Only  a  very  faint  contrast  is  seen  marking  the 
circular  irradiated  area  at  this  fluence.  A  spot  irradiated  at 
higher  fluence  is  seen  in  Figs.  7(c)  and  7(d).  The  faint  re¬ 
gion  here  surrounds  a  dark-ringed  mottled  region.  The  cir¬ 
cular  shapes  of  the  features  are  the  result  of  the  circularly 
symmetric  Gaussian  spatial  profile  of  the  laser  beam.  If  a 
particular  damage  feature  of  radius  r  can  be  associated  with  a 
fluence  threshold /,h ,  this  radius  can  be  related  to  the  aver¬ 
age  pulse  fluence / and  the  Gaussian  radius  <7.  The  average 
pulse  fluence  is  defined  as  the  pulse  energy  E  divided  by  the 
Gaussian  area /=  E  Ar.  The  radius  r  is  thus  given  by 


F - d - H 


FIG  5  Schematic  diagram  of  the  layer  configuration  used  in  the  reflectivity 
calculations  of  partially  molten  films 


FIG.  6  Reflectivity  of  partially  molten  100-nm-thick  films,  calculated  as  a 
function  of  the  thickness  of  the  liquid  layer  (</, )  for  the  geometry  of  Fig  5 
The  angle  of  incidence  is  10"  (a)  Pure  gold  film;  (b)  multilayer  consisting 
of  alternating  10-nm  layers  of  copper  and  gold 


/„  (1) 

or 

r2  =  <r  [  In  (/)  —  In (/,h  )  ]  .  (2) 

The  solid  line  in  Fig.  8  is  a  fit  to  this  function  for  the  sharp 
ring  feature,  one  of  which  is  shown  in  Fig.  7(c).  The  fluence 
threshold  for  this  feature  corresponds  to  the  beginning  of 
region  III  in  Fig.  2,  the  permanent  damage  region.  The  dot¬ 
ted  line  in  Fig.  8  corresponds  to  the  radius  of  the  faint  region 
seen  in  Figs.  7(a)  and  7(c).  This  fluence  threshold  corre¬ 
sponds  to  the  beginning  of  region  II  in  Fig.  2,  the  melting 
regime. 


III.  EXPERIMENTAL  RESULTS 


Figure  2  displays  the  change  in  reflectivity  of  a  copper/ 
gold  multilayer  for  a  particular  pump-probe  delay  at  a  var¬ 
iety  of  fluences.  After  collecting  the  reflectivity  versus 
fluence  data  for  a  particular  time  delay,  the  delay  line  was 
changed  and  the  measurements  were  repeated.  The  reflectiv¬ 
ity  versus  fluence  data  were  then  collected  into  fluence  bins, 
averaged,  and  plotted  as  reflectivity  versus  time  for  a  fluence 
bin  in  the  melting  fluence  regime  II.  Such  a  plot  is  shown  in 
Fig.  9  for  a  30-nm-thick  gold/copper  multilayer  with  10-nm 
repeat  length  on  a  glass  substrate.  The  film  thickness  is  of  the 
order  of  the  optical  absorption  depth.  Thus,  the  drop  in  re- 
P  activity  can  be  associated  with  the  melting  of  the  entire 
film,  and  the  recovery  of  the  reflectivity  back  to  its  original 
value  with  the  movement  of  the  crystal/melt  interface  from 
the  substrate  to  the  free  surface.  The  interface  velocity  u  is 
given  by 


=  d_  =  300x10' 10 
U~  t  4000x  IQ- 12 


=  7.5  m/s , 


(3) 


where  d  is  the  film  thickness  and  t  the  recovery  time.  This  is 
still  a  relatively  low  velocity,  due  to  the  poor  thermal  con- 
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(a) 


(b) 


FIG.  7.  Optical  micrographs  of  laser-irradiated  areas 
for  the  sample  of  Fig  4(a)  Same  fluencc  as  in  Fig.  4; 
bright  field  (b)  Dark  field  of  (a),  (c)  Higher 
fluence;  bright  field  (d)  Dark  field  of  (c). 


(C) 


(d) 


ductivity  of  the  glass.  Note  that  the  maximum  change  in 
reflectivity  in  Fig.  9  is  less  than  3%.  The  fluence  threshold 
for  melting  is  lower  than  for  the  sample  in  Fig.  2  because  of 
the  smaller  film  thickness.  Measurements  in  the  nanosecond 
regime,  such  as  those  by  Tsao  et  al.  on  aluminum,  have 
shown  that  velocities  deduced  from  changes  in  reflectivity 
are  in  agreement  with  those  obtained  from  changes  in  the 
conductivity  in  a  partially  molten  sample.16  The  transient 
conductivity  technique,  however,  cannot  readily  be  ex¬ 
tended  to  the  picosecond  regime. 

The  change  in  reflectivity  as  a  function  of  time  for  a  20- 
nm  gold  film  on  a  sapphire  substrate  is  shown  in  Fig.  10.  The 
reflectivity  recovery  began  at  400  ps  and  was  completed  by 
600  ps.  Equation  ( 3 )  yields  a  velocity  of  1 00  m/s  in  this  case. 
This  is  the  highest  velocity  for  growth  of  a  crystal  from  its 
melt  ever  measured  by  any  technique.  This  value  is  probably 


FIG  8  Square  of  the  radius  of  the  damage  features  of  Fig.  7(c )  as  a  function 
of  the  logarithm  of  the  incident  laser  fluence.  The  solid  line  is  a  a  least- 
squares  fit  for  the  sharp  ring  The  dotted  line  is  a  fit  with  the  same  slope  for 
the  radius  of  the  faint  white  region 


still  only  a  lower  limit.  The  value  of  the  regrowth  time  used 
here  is  the  difference  between  the  last  time  at  which  the  re¬ 
flectivity  is  low  and  the  first  at  which  it  has  relumed  to  its 
initial  value;  the  actual  recovery  time  may  be  smaller.  Fur¬ 
thermore,  if  the  interface  is  not  fun.,  planar,  the  recovery  of 
the  reflectivity  would  be  smeared  out  in  time  so  that  the 
actual  velocity  would,  again,  be  larger  than  the  value  deter¬ 
mined  above.  On  the  other  hand,  if  the  melt  front  did  not 
reach  the  back  surface,  the  above  value  would  be  too  large.  If 
this  were  the  case,  however,  the  reflectivity  trace  of  Fig.  10 
should  have  a  V  shape;  instead,  it  has  a  flat  portion  (between 
200  and  400  ps)  that  corresponds  to  the  time  that  the  entire 
film  was  fully  melted. 

The  time-resolved  reflectivity  technique  was  also  ap¬ 
plied  to  a  bulk  sample.  The  resulting  data  showed  more  scat¬ 
ter,  however,  due  to  spot-to-spot  variation  of  its  surface 
characteristics.  A  plot  of  reflectivity  versus  time  is  shown  in 
Fig.  11  for  a  bulk  copper  sample.  In  this  case  it  is  necessary  to 


FIG  9  Relative  reflectivity  change  as  a  function  of  time  for  a  30-nm-thick 
gold/copper  multilayer  with  10-nm  repeat  length  on  a  glass  substrate,  irra¬ 
diated  with  a  4-mJ/cm’  laser  fluence 
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FIG.  10.  Relative  reflectivity  change  as  a  function  of  time  for  a  20-nm-thick  FIG.  12  Relative  changes  in  reflectivity  and  transmittivity  as  a  function  of 

gold  film  ona  sapphire  substrate,  irradiated  with  a  24-mJ/cm’ laser  fluence.  llme  f°r  an  aluminum  film  on  a  sapphire  substrate  (/?„  =  0  6,  T„  ~  0  2), 

irradiated  with  a  23-mJ/cnv  laser  fiuence. 


associate  the  reflectivity  recovery  time  with  the  movement  of 
the  interface  over  the  optical  absorption  depth,  25  nm  in  this 
material.  This  gives  a  velocity  of  60  m/s  for  the  crystal  re¬ 
growth.  It  might  be  expected  that  the  regrowth  velocity  for  a 
bulk  metal  would  be  higher  than  that  for  a  metal  film  on 
sapphire  because  of  the  higher  thermal  conductivity  of  the 
metal  substrate.  However,  the  melt  in  the  bulk  experiment 
was  probably  somewhat  deeper,  and  the  crystallization  front 
may  have  slowed  due  to  latent  heat  rejection  by  the  time  it 
reached  the  25-nm  surface  region  where  crystallization 
could  be  detected. 

In  all  cases,  the  decrease  in  reflectivity  associated  with 
melting  was  completed  in  less  than  100-200  ps.  This  is  also 
shown  for  the  particular  case  of  an  aluminum  film  on  a  sap¬ 
phire  substrate  in  Fig.  12.  These  melting  times  are  in  appar¬ 
ent  disagreement  with  the  picosecond  electron  diffraction 
experiments  of  Williamson,  Mourou,  and  Lin  on  free-stand¬ 
ing  20-nm-thick  aluminum  films,  in  which,  at  low  fluence  ( 8 
mJ/cm2  absorbed),  no  melting  was  observed  until  1  ns  after 
irradiation.  Rapid  melting,  similar  to  that  observed  here, 
was  also  seen  in  their  experiments,  but  at  higher  fluences  ( 1 3 
mJ/cm2  absorbed).  It  can  be  argued"1  that  if  melting  re¬ 
quires  1  ns  to  be  observed  at  the  lower  fluences,  it  would  not 
occur  at  all  in  the  experiment  described  in  this  paper,  since 
by  that  time  all  the  deposited  heat  would  have  drained  away 
through  the  substrate  and  cooled  the  surface  to  below  the 
melting  temperature.  Thus,  the  difference  in  sample  geome¬ 
try  may  have  led  to  an  apparent  higher  threshold  fluence  for 
melting  in  the  present  experiments.  The  actual  absorbed 
fluence  for  the  sample  of  Fig.  10  is  about  5  mJ/cm2,  at  the 


TIME  (pj) 


FIG.  II.  Relative  reflectivity  change  as  a  function  of  time  for  a  bulk  copper 
sample,  irradiated  with  a  ZJ-mJ/cm’  laser  fluence. 


low  end  of  the  values  used  in  the  experiments  of  Williamson 
and  co-workers.  The  film  thickness  for  our  sample,  however, 
is  probably  smaller  than  20  nm,  and  no  doubt  has  somewhat 
different  optical  properties,  all  of  which  makes  a  direct  com¬ 
parison  of  the  fluences  in  the  two  experiments  difficult. 

IV.  DISCUSSION  AND  COMPUTER  MODELING 

An  analytic  calculation  of  the  temperature  at  each  point 
in  the  sample  under  the  conditions  of  pulsed  laser  quenching 
is  very  complex,  especially  if  the  temperature  dependence  of 
the  material  properties,  the  kinetics  of  the  various  phase 
transitions,  and  the  mechanical  effects  must  be  taken  into 
account.  In  order  to  make  the  problem  tractable,  it  is  neces¬ 
sary  to  make  a  series  of  rather  restrictive  approximations. 
The  main  assumption  is  that  the  Laplace  equation  is  a  valid 
description  of  the  heat  flow  process,  even  in  the  presence  of 
large  thermal  gradients, 

4r  =  i-V(/rVD+  -ie(f,t),  (4) 

dt  C  C  C 

where  T  is  the  temperature,  C  the  heat  capacity  per  unit 
volume,/  the  laser  heat  source,  r  the  position  in  the  film,/  the 
time,  a  the  thermal  conductivity,  and  Q  the  heat  loss  to  var¬ 
ious  sinks  as  well  as  the  latent  heat  released.  One  of  the  as¬ 
sumptions  implicit  in  the  use  of  the  Laplace  equation  is  that 
convection  either  does  not  occur  or  is  negligible.  Calcula¬ 
tions  indicate  that  this  is  reasonable  because  of  the  small 
thickness  of  the  melt  and  its  short  duration.19  It  has  been 
demonstrated20  that  a  one-dimensional  equation  adequately 
describes  the  case  of  a  circularly  symmetric  Gaussian  beam 
on  a  laterally  homogeneous  material  if  the  beam  radius  is 
greater  than  twice  the  thermal  diffusion  length,  a  condition 
that  is  well  satisfied  in  the  present  experiments.  Possible  heat 
losses  include  conduction  to  the  sample  holder  or  the  sur¬ 
rounding  air,  blackbody  radiation,  and  evaporation  of  the 
material  from  the  surface.  These  two  types  of  conduction  are 
negligible  for  the  time  of  the  experiment.  Blackbody  radi¬ 
ation  has  been  shown  to  be  negligible  as  well.21  Evaporation 
losses  were  included  in  the  computer  model  but  were  also 
found  to  be  very  small,  in  agreement  with  earlier  simple  esti¬ 
mates.2223 

Most  of  the  existing  computer  simulations  of  pulsed  la¬ 
ser-induced  transformations  are  based  on  the  model  by 
Baeri.24  In  this  model  the  differential  equation  describing  the 
heat  flow  is  replaced  with  a  set  of  finite  element  equations. 


133 


J.  Appl  Phys  .  Vol  65.  No  1.1  January  1989 


MacHrinald  Maivo?7*  anti  c;n:*ooor> 


The  sample  is  divided  into  slices  of  width  A z,  and  time  into 
intervals  A/.  The  material  parameters  can  be  made  into  func¬ 
tions  of  both  the  temperature  and  the  phase  of  the  material. 
When  the  temperature  of  the  /th  layer  reaches  the  melting 
point,  the  temperature  change  A  Tis  held  equal  to  zero  until 
the  accumulated  heat  input  is  equal  to  the  latent  heat  of 
fusion.  The  intervals  Azand  A/arechosen  tofulfil!  the  stabil¬ 
ity  requirements  for  numerical  convergence:  kA/  <  C(  Az)  V 
2.  This  model  has  proved  to  be  in  good  agreement  with  a 
variety  of  experimental  results.  25'2* 

In  this  work  a  similar  technique  was  employed,  based  on 
theCrank-Nicholson  method  for  solving  the  heat  flow  equa¬ 
tion.29  In  this  alogrithm  a  matrix  equation  is  set  up  that 
relates  the  temperature  in  each  layer  to  the  value  in  the 
neighboring  layers,  and  this  entire  matrix  equation  is  solved 
for  each  new  time  interval.  Such  an  implicit  algorithm  has 
the  advantage  that  there  are  no  numerical  stability  restric¬ 
tions  on  the  size  of  the  time  interval  for  a  given  distance 
interval.  Approximations  must  be  made  for  the  coefficients, 
since  they  are  a  function  of  the  material  parameters  at  the 
unknown  temperature.  A  two-step  process'0  was  employed. 
Material  parameters  determined  from  the  previous  tempera¬ 
tures  were  used  in  the  first  pass  to  calculate  new  estimated 
temperatures.  These  estimated  temperatures  were  used  to 
compute  new  material  parameters,  which  were  then  used  in 
the  second  pass  to  calculate  better  values  for  the  new  tem¬ 
peratures.  This  procedure,  which  must  include  an  estimate 
of  the  latent  heat  released  by  the  moving  interface,  intro¬ 
duces  a  new  restriction  on  the  maximum  size  of  the  time  step 
for  numerical  stability.  However,  it  is  linear  in  Az  and  there¬ 
fore  less  severe  than  in  the  previous  model. 

For  the  case  of  a  one-dimensional  homogeneous  sample 
of  infinite  thickness,  with  temperature-independent  material 
parameters,  irradiated  at  constant  power  for  finite  duration, 
the  equation  can  be  solved  analytically.31  The  computer 
model  solution  for  this  case  agrees  well  with  the  exact  solu¬ 
tion. 

The  most  important  difference  between  the  present 
model  and  its  predecessors  is  the  incorporation  of  the  actual 
kinetics  of  melting  and  crystallization.  Unlike  in  the  nano¬ 
second  or  longer  pulse  experiments,  for  which  most  of  the 
earlier  models  were  developed,  in  picosecond  experiments 
the  interface  temperature  during  regrowth  can  be  several 
hundred  degrees  below  the  melting  temperature,  and  the  ve¬ 
locity  of  the  interface  can  therefore  be  significantly  slower 
than  that  of  the  melting  temperature  isotherm. 

It  is  clear  that  the  kinetics  of  crystal  growth  depend  on 
the  nature  of  the  solid-liquid  interface  and  the  details  of  the 
atomic  movements,  as  well  as  on  the  macroscopic  heat  flow. 
The  crystallization  rate  can  be  described  as  the  difference 
between  a  freezing  and  a  melting  flux.  In  equilibrium,  the 
two  fluxes  are  equal.  The  fluxes  can  be  evaluated  using  tran¬ 
sition  state  rate  theory.'2"  It  is  assumed  that  there  exists  an 
activated  state  with  a  free  energy  AG  *  above  that  of  the  li¬ 
quid.  If  v  is  the  atomic  vibration,  or  attempt  frequency,  the 
number  of  activated  configurations  per  interfacial  atom 
formed  per  unit  time  from  the  liquid  or  the  solid  is,  respec¬ 
tively,  ve  -  *°’/KT  or  ve  ~  liC’  ■*  ac,/*r,  where  AG  is  the  free 
energy  of  crystallization.  The  growth  velocity  u  is  then  ob¬ 


tained  from  the  difference  between  these  two  rates 

u=/vMe  *cv*r-e  -|4C"4G,/*r)  .  (5) 

where  A  is  the  average  jump  distance  (on  the  order  of  an 
interatomic  distance )  and /  the  fraction  of  interfacial  sites  at 
which  attachment  can  occur.  This  fraction  is  near  unity  for 
“rough”  crystal-melt  interfaces,  such  as  those  in  metals,  but 
is  small  for  “smooth”  interfaces,  such  as  the  low-index  ones 
in  silicon  at  small  undercooling.3 

It  is  useful  to  consider  the  atomic  mechanisms  associat¬ 
ed  with  particular  values  of  the  free  energy  of  activation 
AG*.  If  AG*  is  nonzero,  the  growth  velocity  initially  in¬ 
creases  with  undercooling  due  to  the  increase  in  driving  free 
energy  AG,  but  eventually  it  falls  off  due  to  the  decrease  in 
atomic  mobility.  A  particular  example  is  diffusion-limited 
growth,  in  which  the  motion  of  an  atom  from  the  liquid  to  the 
crystal  requires  a  rearrangement  similar  to  that  for  a  diffu¬ 
sive  jump  in  the  liquid,  so  that  AG  *  then  corresponds  to  the 
free  energy  of  activation  for  the  liquid  diffusivity.  This  type 
of  growth  occurs,  for  example,  in  concentrated  metallic  al¬ 
loys  or  in  systems  with  a  large  crystalline  unit  cell.  If  the 
interfacial  temperature  falls  below  the  diffusional  freezing 
temperature  7j, ,  crystal  growth  stops  and  the  remaining  liq¬ 
uid  becomes  a  glass.  It  should  be  kept  in  mind,  however,  that 
some  very  dilute  amorphous  metallic  alloys  produced  by  va¬ 
por  deposition  on  a  cold  substrate  can  crystallize  at  tempera¬ 
tures  as  low  as  30  K.  u  This  is  well  below  Tg ,  indicating  that 
crystallization  in  such  systems  proceeds  by  some  faster,  non- 
diffusional,  process,  such  as  the  collision-limited  one  dis¬ 
cussed  below.  The  maximum  growth  velocity  in  the  diffu¬ 
sion-limited  regime  uD  can  be  estimated  from  the  liquid 
diffusivity,  which  in  metals  is  DL  ~  10~ 5  cm2/s, 

uD-~A /(A2/DL)  =  (DL/A)  —  3  m/s,  ^  (6) 

since  {A  2/DL )  is  the  time  required  to  diffuse  an  interatomic 
distance  A. 

If  AG  *  is  zero,  the  frequency  with  which  atoms  cross  the 
interface  is  the  thermal  vibration  frequency.  This  corre¬ 
sponds  to  collision-limited  growth,  in  which  the  maximum 
velocity  is  the  speed  of  sound.  It  has  been  proposed  that  this 
mechanism  governs  the  crystallization  of  pure  metals.35  36 
Perhaps  the  best  evidence  for  this  hypothesis  is  Walker’s 
measurements"  of  the  velocity  of  nickel  dendrites  grow  ing 
into  an  undercooled  melt  at  apparent  velocities  up  to  50  m/s. 
An  analysis  by  Coriel)  and  Turnbull'*  indicates  that  the  data 
is  best  fit  by  a  model  that  incorporates  collision-limited 
growth.  There  is  also  indirect  evidence  from  a  laser  quench¬ 
ing  experiment.  Lin  and  Spaepen'9  observed  that  a  sputter- 
deposited  iron  sample  with  an  initial  grain  size  of  50  nm 
contained  grains  with  a  diameter  as  large  as  1  // m  after  irra¬ 
diation  with  a  30-ps  pulse.  Using  a  simple  heat  flow  estimate 
for  melt  duration  of  —  2  ns  yields  an  apparent  lateral  growth 
velocity  of  several  hundred  meters  per  second,  which  is 
much  faster  than  the  maximum  diffusion-limited  growth  ve¬ 
locity. 

Since  the  velocities  measured  in  the  present  experiments 
are  also  much  greater  than  the  maximum  velocity  in  diffu¬ 
sion-limited  growth,  the  model  calculation  is  done  with  col¬ 
lision-limited  growth  kinetics.  Furthermore,  since  it  is  ap- 
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plied  to  metallic  systems,  the  site  factor  /is  set  equal  to  unity. 
With  the  approximation  AG  =  AS(T„  —  T),  where  AS  is 
the  entropy  of  fusion  and  T„  is  the  melting  point,  Eq.  (5) 
becomes 


u  =  u,  ( 1  —  e 


tss/R) j  i  t„  -  n/r] 


)  , 


(7) 


where  us  is  the  speed  of  sound.  For  the  case  of  small  devia¬ 
tions  from  equilibrium  this  becomes 

u  =  us(AS/R7~)(  Tm  -  T)  ,  (8) 


yielding  the  familiar  result  that  the  growth  velocity  is  pro¬ 
portional  to  the  undercooling  near  equilibrium. 

It  should  be  kept  in  mind  that  the  transition  state  rate 
theory  is  a  phenomenological  one,  which  cannot  be  expected 
to  apply  rigorously  under  all  conditions.  In  particular,  the 
extension  of  the  theory  to  the  condition  AG*— 0  (t.e.,  the 
removal  of  the  activation  barrier)  is  open  to  question. 
Further,  it  can  be  seen  that  the  absolute  value  of  the  melting 
velocity  in  Eq.  (7)  is  larger  than  the  speed  of  sound  as  the 
temperature  goes  to  positive  infinity.  This  seems  physically 
implausible.  Tsao elal.4"  resolve  this  problem  by  considering 
separately  the  entropic  and  enthalpic  contributions  to  the 
free  energy  barrier.  By  taking  the  entropy  as  the  barrier  to 
crystallization  and  the  enthalpy  as  the  barrier  to  melting, 
they  obtain 

u  =  ut(e~  -  e~ ,A"  *n) ,  (9) 


where  AH  is  the  enthalpy  of  fusion.  Since  AH  —  T^AS,  this 
can  be  w  ritten  as 


u  -  ute 


- (US/Rl 


(l-e 


is s -Rt|  (  rv  -  n/r] 


)  , 


(10) 


in  which  the  speed  of  sound  is  now  the  limit  for  the  melting 
velocity,  and  the  speed  of  sound  reduced  by  the  factor 
exp(  —  AS /R)  for  the  crystallization  velocity.  Linearizing 
this  equation  yields 


u  =  u, 


AS  ( T -  T„)  e_AS/M 

R  T„ 


This  is  the  equation  used  in  the  computer  model. 

The  mode!  was  tested  by  applying  it  to  nanosecond 
pulses  on  silicon,  and  comparing  the  results  to  those  ob¬ 
tained  with  Thompson’s  computer  code,27  which  was  based 
on  Baeri’s  model  described  above,  with  good  agreement. 
This  version  of  Thompson’s  code,  which  requires  a  smaller 
time  step  and  does  not  include  an  interface  response  func¬ 
tion,  is  not  expected  to  be  accurate  in  the  picosecond  regime, 
but  has  been  found  to  agree  well  with  nanosecond  results  on 
silicon.27 

To  model  the  present  experimental  results,  a  computa¬ 
tion  was  made  for  a  20-nm  gold  film  on  a  sapphire  substrate 
irradiated  with  a  20-ps  pulse,  using  the  materials  param¬ 
eters41  listed  in  Table  I.  The  melt  depth  and  the  correspond¬ 
ing  interface  velocity  were  calculated  as  a  function  of  time, 
and  are  shown  in  Fig.  13.  The  model  predicts  velocities  on 
the  order  of  100-200  m/s,  for  a  fiuence  very  near  the  melting 
threshold.  A  larger  fiuence  gives  a  longer  melt  duration  and 
a  lower  velocity.  Figure  13  is  in  general  agreement  with  the 
experimental  results  in  Fig.  10.  An  exact  comparison  is  diffi¬ 
cult  because  of  the  uncertainty  in  many  of  the  materials  pa¬ 
rameters  at  high  temperatures,  particularly  in  the  absorbti- 
vity  and  reflectivity  at  the  pump  pulse  wavelength. 
Calculations  for  a  bulk  copper  specimen  gave  similar  results. 


V.  CONCLUSIONS 

The  experiments  described  in  this  paper  demonstrated 
that  pump-probe  reflectivity  measurements  can  be  used  to 
make  time-resolved  studies  of  melting  and  crystallization  in¬ 
duced  by  picosecond  pulsed  laser  irradiation  of  simple  met- 


TABLE  I  Parameters  used  in  the  simulations 


Comparison 
with  analytic 
solution 

Silicon 

Copper 

Gold 

Sapphire 

«t,  solid  (J/cm  s  K) 

40 

1  68 

40 

3  18 

0.3 

C,  solid  (J/cm'  K)  —  A  ■¥  BT  +  D/T 

A 

3  04 

1.17 

3  19 

2  52 

4.25 

B 

0  209 

-  0  885 

0.072 

0.714 

D 

-  0  217 

a,  solid  (cm1) 

ICT 

1 340  Xesp(  77427) 

4  84  x  10' 

4  84 v 10' 

001 

R,  solid 

06 

0.337 

06 

06 

F  (J/cm!) 

0  5 

10 

0.07 

0.06 

7.,,.r  (J/mol) 

4  3x  I04 

14  5159 

4  3  «  10' 

3.2  :  10' 

T.  (K) 

2867 

2  673 

2867 

3080 

p  (g/cm1) 

8  9b 

2  33 

8  96 

19.32 

3  97 

M  (g/mol) 

63.5 

63  5 

28.1 

197 

102 

T„  (K) 

1  687 

1356 

1336 

2300 

L  (J/cm') 

2  270 

1836 

v  1243 

4230 

A S/R 

3  6 

1  12 

111 

5.33 

C,  liquid  (J/cm'  K ) 

2  56 

4  43 

2.87 

k,  liquid  (J/cm  s  K ) 

1  4 

40 

3.18 

R.  liquid 

0  72 

0  5 

0  5 

a.  liquid  (cm1) 

69  325 

4  8  >10' 

4  8  >10' 
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FIG.  13.  Model  calculations  of  the  melt  depth  and  interface  velocity  as  a 
function  of  time  for  a  20-nm  gold  film  on  a  sapphire  substrate.  The  simula¬ 
tion  parameters  are  listed  in  Table  1. 

als,  such  as  gold,  copper,  or  aluminum,  even  though  the  dif¬ 
ference  between  their  reflectivities  in  the  liquid  and 
crystalline  is  small. 

That  the  observed  reflectivity  changes  corresponded  to 
melting  was  demonstrated  by  performing  the  experiments 
on  artificial  multilayers  of  fully  miscible  metals,  such  as  cop¬ 
per  and  gold,  and  confirming  that  melting  occurred  by  ob¬ 
serving  the  leveling  of  the  composition  profile  as  a  result  of 
liquid  phase  diffusion  by  Auger  spectroscopy. 

Crystallization  velocities  up  to  100  m/s  were  measured. 
These  are  the  largest  directly  measured  values  for  any  mate¬ 
rial.  It  is  clear  that  the  atomic  scale  mechanism  governing 
crystallization  in  these  simple  metals  cannot  be  purely  diffu- 
sional  in  nature,  since  the  measured  velocities  are  substan¬ 
tially  greater  than  the  maximum  growth  velocity  in  the  diffu¬ 
sion-limited  growth  regime.  On  the  other  hand,  these 
experiments  lend  strong  support  to  the  idea  that  crystalliza¬ 
tion  in  these  systems  is  collision  limited,  since  they  could  be 
accounted  for  satisfactorily  by  a  computer  model  that  incor¬ 
porated  both  the  heat  flow  processes  a.,d  collision-limited 
kinetics. 
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1.  Introduction 

The  rapid  extraction  of  heat  from  a  metallic  melt  can  lead  to  a  substantial 
undercooling  below  the  thermodynamic  crystallization  temperature .  In  this  paper 
some  of  the  fundamental  kinetic  processes  of  crystallization  at  large  undercooling 
are  reviewed:  homogeneous  crystal  nucleation,  high-velocity  crystal  growth  induced 
by  pulsed  laser  irradiation,  and  the  formation  of  metastable  crystalline  phases.  An 
extended  version  of  this  paper  can  be  found  in  two  recent  publications  by  the 
authors  [1,2]  . 

2.  Homogeneous  Nucleation 

a .  Experimental  observations. 

Using  dispersion  and  fluxing  techniques  Turnbull  and  his  associates  demonstrated  in  the 
early  1950's  that  AT0fr  =  AT0/Tj  (AT0:  undercooling  at  measurable  homogeneous 

nucleation;  Tj:  liquidus  temperature)  was  at  least  -0.18  for  most  of  the  close-packed 

and  b.c.c.  metals  [3]  and  -0.32  for  Hg  [4].  At  the  time  the  large  magnitudes  of  these 
undercoolings  were  quite  unexpected  because  of  the  then  prevalent  view  that  the  atomic 
short  range  order  (SRO)  in  liquid  metals  was  essentially  similar  to  that  in  close- 
packed  crystals.  Less  surprising  were  the  results  of  Vonnegut  [5],  Turnbull,  and 
others  that  the  AT0,  r  in  liquids  that  crystallize  to  structures  with  a  coordination 
number  less  than  eight,  also  were  0.15-0.20,  at  least. 

Later  investigators  using  similar  techniques  found  onset  undercoolings  greater  than  0.3 
for  some  of  these  liquids  and  between  0.2  and  0.3  for  several  close  packed  metals. 

More  particularly,  Bosio  et  al.  [6]  reported  AT0>r  approaching  0.5  for  certain 
emulsions  of  Ga  droplets.  Perepezko  and  coworkers  [7]  obtained  AT0>  r  in  the  range  of 

0.3  to  0.4  for  droplet  emulsions  cf  In,  Sn  and  Bi .  Recently  Devaud  and  Turnbull  [8] 
observed  that  some  0.3-0. 5  mm  diameter  Ge  droplets,  fluxed  in  molten  dehydrated  B2O3, 

exhibited  ATQ, r  *  0.35.  Ovsienko  [9],  Gomersall  et  al.  [10]  reported  AT0( r  ranging 
between  0.23  and  0.27  for  droplets  of  Ge,  Co,  Ni .  Table  I  summarizes  some  of  the  more 
recent  data  on  onset  undercooling. 

It  is  generally  presumed  that  the  undercoolings  at  the  onset  of  measurable  homogeneous 
nucleation  must  equal  or  exceed  the  maximum  reported  experimental  values,  properly 


adjusted  for  cooling  rate  and  specimen  size.  However,  there  are  conditions,  as 
explained  by  Turnbull  [12],  under  which  this  presumption  may  not  be  correct.  First,  it 
is  evident  that,  when  comparing  results  of  different  investigators,  proper  allowance 
must  be  made  for  experimental  uncertainties  in  the  calibration  of  TQ,  the  temperature 

at  nucleation  onset,  relative  to  the  liquidus,  Tj.  However,  apart  from  calibrational 
disparities  there  are  other  reasons  why  there  might  be  some  apparent  or  actual 
specimen-to-specimen  dispersion  of  AT0>r  values  at  homogeneous  nucleation. 


Table  I.  Reported  scaled  undercoolings,  AT0/r 
(except  for  Hg  uncorrected  for  differences  in  cooling  rate) 
at  measured  nucleation  onset  for  several  pure  metals. 


Metals 

AT0,r 

Specimen  Diameter 

Ref. 

Fe 

0.23  ±  0.01 

~  0.7  cm 

10 

Co 

0.27  ±  0.01 

60-200  Jim 

9 

Ni 

0.27  ±  0.01 

80-200  pm 

9 

In 

0.26 

2-5  pm 

7 

Hg 

0.33  -  0.35 

2-5  pm 

4,11 

Ga 

0.50  ±  0.02 

2-5  pm 

6 

Ge 

0.35  ±  0.02 

0.05  cm 

8 

Sn 

0.37 

2-5  pm 

7 

Bi 

0.42 

2-5  pm 

7 

An  apparent  dispersion  may  arise  between  various  emulsions  of  droplets  coated  by  solid 
films,  because  substantial  hydrostatic  stresses  can  develop  in  the  smaller  droplets 
(diamater  <2-3  microns)  as  they  are  undercooled,  due  to  the  differential  volume 
response  of  the  liquid  relative  to  its  coating  film.  Consequently,  the  apparent  ATq> r 
could  be  displaced  considerably  from  that  which  would  be  exhibited  by  the  unstressed 
liquid.  An  analysis  by  the  authors  [13]  indicated  that  this  displacement  could  be  as 
large  as  +10%  for  small  Afi  droplets  coated  by  thin  films  of  their  oxide.  In  view  of 
the  high  mechanical  strength  of  solids  in  thin  film  form  and  the  great  resistance  of 
liquid  metals  to  cavitation,  the  possible  importance  of  such  stress  effects  in 
experiments  on  small  droplets  cannot  be  lightly  dismissed. 

An  actual  dispersion  of  ATQf  r  for  homogeneous  nucleation  can  result  from  specimen-to- 
specimen  variation  in  the  concentrations  of  those  homophase  impurities  which  have  a 
lower  standard  free  energy  of  solution  in  the  crystalline  than  in  the  liquid  phase  of 
the  metals  [12]  .  Such  impurities,  even  in  amounts  of  order  of  a  few  at  ppm,  can, 
without  measurably  displacing  Tj,  reduce  AT0>r  considerably  from  its  pure  metal  value. 

This  "chemical"  effect  on  nucleation  is  distinct  from  that  of  the  heterophase 
impurities  which  lead  to  heterogeneous  nucleation. 

At  steady  state,  isothermal  homogeneous  nucleation  should  be  stochastic  in  volume  and 
time.  This  condition  requires  that  in  droplet  emulsions  the  nucleation  probability, 
k^,  in  a  droplet  of  diameter,  d,  be  proportional  to  Id^,  where  I  is  the  nucleation 
frequency /volume-time .  If  the  undercooling  is  chosen  so  that  the  time  for  crystal 
growth  is  negligible  compared  with  that  for  nucleation  k^  becomes,  virtually,  the 
crystallization  frequency.  Then  the  isothermal  kinetics  of  crystallization  of  droplet 
emulsions  can  be  analyzed  to  ascertain  the  dependence  of  k^  on  d.  In  one  method  (a) 
developed  by  Turnbull  [4]  and  described  in  detail  by  Pound  [14],  the  size  dispersion  of 
the  droplet  population  is  exploited  to  discriminate  between  the  possibilities. 


However,  a  more  conclusive  discrimination  is  achieved  by  method  (b)  in  which  k<j,  where 
d  is  the  diameter  of  the  droplet  of  average  volume,  are  compared  for  separate  emulsions 
with  substantially  differing  d;  then  k^/k^  “  <  /  <^2 )  ^  if  nucleation  is  homogeneous. 

Turnbull's  measurements  of  the  isothermal  rates  of  solidification  of  Hg-Laurate  (Hg-L) 
coated  Hg  droplet  dispersions  [4],  as  well  as  those  with  Cormia  [15]  on  n-hepta-  and  n- 
octadecane  dispersions,  fulfilled  both  tests  (a)  and  (b)  for  volume  nucleation.  Thus 
at  a  given  undercooling  the  rates  were  satisfactorily  described  by  a  single  nucleation 
frequency /volume .  The  onset  undercooling  of  Hg  in  these  experiments  was  -78°  (ATC  r  * 
0.33).  Somewhat  higher  onset  undercoolings  in  differently  prepared  Hg  dispersions  were 
reported  by  Monge  [16]  and  by  Perepezko  and  Rasmussen  [11]  .  When  adjusted  to  the 
cooling  rate  of  Turnbull's  experiments  the  Monge  result  is  about  3°  and  that  of 
Perepezko  and  Rasmussen  about  4-5°  greater  than  obtained  by  Turnbull.  These 
differences  are  not  large  and  could  reflect  only  the  temperature  calibration 
uncertainties  in  the  separate  studies.  Alternatively,  they  could  be  due  to  small 
"chemical"  effects  [12]  as  discussed  in  the  foregoing.  Thus  the  evidence  that 
homogeneous  nucleation  was  exhibited  in  Turnbull's  experiments  on  HgL  coated  Hg 
droplets  seems  quite  conclusive. 

The  temperature  dependence  of  I  in  Turnbull's  experiments  on  Hg  was  satisfactorily 
described  by  the  equation  derived  from  classical  nucleation  theory  but  with  a  prefactor 
about  107  times  too  large  [4],  when  a  and  the  entropy  of  crystallization/volume,  ASV, 

were  assumed  constant.  It  was  shown  elsewhere  [12]  that  this  disparity  could  result 
from  the  chemical  effect  noted  above.  Alternatively,  it  might,  as  Turnbull  suggested 
earlier  [4],  arise  from  temperature  variations  of  ASV  and  a.  The  temperature  variation 
of  ASV(T)  in  the  undercooled  regime  can  be  estimated  roughly  from  heat  capacity,  Cp(T), 
measurements  [1] .  Putting  ASV(T)  into  the  classical  expression  for  I,  with  the 
theoretical  prefactor,  the  Hg  data  can  then  be  fit  within  experimental  error,  with  an 
interfacial  tension  increasing  by  0.10  ergs/cm2-deg  with  temperature.  Over  the 
temperature  regime  of  the  I  measurements  the  fitting  <J  could  be  described  by: 

0  ■  11.9  +  0.10T  ergs/cnP  (1) 

where  T  is  the  absolute  temperature.  This  result  is  in  rough  agreement  with  that  of  a 
calculation  of  Mijazawa  and  Pound  [17]  assuming  ACp  -  0. 

Mijazawa  and  Pound  [17]  measured  the  isothermal  rates  of  solidification  of  two 
dispersions  of  Ga  droplets  which  exhibited  onset  undercoolings  -100°,  well  below  the 
maximum  -150°  claimed  by  Bosio  et  al.  [6].  Their  data  on  the  separate  dispersions 
seemed  to  meet  test  (a)  for  volume  nucleation.  However,  they  appear  to  fail  to  meet 
test  (b) .  Thus,  at  AT  «  100°  the  solidification  time  constant  (l/k^p  for  their  first 

dispersion,  in  which  d  «  5  urn,  was  only  3.5  times  greater  than  that  for  their  second 
dispersion  in  which  d  (=  21  [im)  was  about  four  times  larger.  Thus  their  k^“  appeared  to 

be  proportional  to  d,  rather  than  to  dP,  which  would  have  been  required  by  test  (b)  for 
volume  nucleation.  The  origin  of  such  an  unusual  dependence  of  k^“  on  d  is  not  clear. 

b.  Theoretical  developments 

The  crystal-melt  inter facial  tension.  Macroscopically  the  nucleation  resistance 
indicates,  in  terms  of  simple  nucleation  theory,  a  scaled  crystal-melt  interfacial 
tension  at  the  temperature  of  measurable  nucleation  onset  of  about  0.55  for  the  close- 
packed  metals  and  0.60  for  Hg.  This  scaled  tension  is  defined  as: 

a  ■*  v^ f  ^d/Ahf  (2) 


where  v  is  the  volume/atom  and  Ahf  is  the  heat  of  fusion/atom. 

In  metals  the  volume  changes  on  solidification  are  small  and  at  constant  density  the 
energy  is  quite  insensitive  to  the  prexie  atomic  siting.  From  this  it  appears  that  a 
crystal-melt  interface,  without  a  density  deficit  should  be  formable  at  little  energy 
cost.  However,  depending  on  liquid  structure,  formation  of  such  an  interface  may 
require  localization,  with  consequent  entropy  loss,  of  atomic  positions  in  the  liquid 
adjacent  to  the  interface  [18]  . 

Spaepen  and  Meyer  [19,20]  calculated  that  the  contribution  to  the  scaled  tension,  a, 
associated  with  this  localization  at  the  (111)  face  of  an  f.c.c.  crystal  is  -0.86  at 
Tm.  This  value  of  a  should  be  an  upper  limit  and  it  is,  indeed,  well  above  those 
indicated  by  the  crystal  nucleation  kinetics.  Either  the  negentropy  or  the  excess 
energy  of  the  interface  could  be  made  to  vanish  by  particular  choices  of  the 
interfacial  configuration  but  they  cannot  be  made  to  vanish  simultaneously  [18]  .  The 
nature  of  the  interfacial  tension  (negentropic  or  energetic)  does  make  a  difference  for 
the  temperature  dependence.  Specifically,  if  the  interfacial  tension  is  negentropic  in 
nature,  it  is  expected  to  have  a  positive  temperature  coefficient.  This  is  what  is 
observed  in  the  Hg  experiments,  as  shown  by  equation  (1),  and  the  results  can  be 
accounted  for  quantitiatively  by  a  fully  negentropic  model  [21] . 

Microscopically  the  high  onucleation  resistance  of  liquid  metals  implies  that  their 
crystallization  must  be  attended  by  some  essential  reconstruction  of  the  atomic  short 
range  order.  Frank  [22]  suggested  that  this  reconstruction  might  be  from  a 
predominantly  icosahedral  short  range  order  in  the  liquid  to  that  in  the  close  packed 
crystal.  In  support  of  this  suggestion  he  showed  that  if  the  interatomic  interaction 
potential  is  of  the  Lennard-Jones  type,  the  energy  of  a  13  atom  cluster  is  actually 
lower  in  the  icosahedral  than  in  either  of  the  close-packed  crystalline  configurations. 
The  icosahedral  configuration  is  a  polytetrahedral  one  and  includes  no  half-octahedra, 
which  are  essential  components  of  close-packed  crystalline  structures.  It  has  also 
been  established  that  the  energetic  preference  for  polytetrahedral  over  close-packed 
crystalline  configurations  extends  well  beyond  13  to  several  tens  of  atoms  clusters  and 
that  the  density  in  these  polytetrahedral  clusters  is  higher  than  in  the  crystalline 
ones  [23]  .  In  the  dense  random  packed  structure  of  uniform  hard  spheres  it  is  found 
that  tetrahedral  configurations  predominate  while  there  is  a  conspicuous  paucity, 
relative  to  the  crystals,  of  half-octahedral  ones  [24] .  The  appearance  of  a 
polytetrahedral  type  structure  in  random  packing  suggests  that  it  is  entropically 
favored  and  that  it  may  be  the  most  likely  prototypical  model  structure  for  monatomic 
liquids.  The  entropy  of  a  condensed  system  should  increase  with  the  range  of  local 
density  fluctuations  permissible  at  a  given  energy.  Because  of  the  high  density  of 
polytetrahedral  clust  •••'-s  this  range  should  be  substantially  greater  in  random  packed 
structures  than  in  those  where  the  crystalline  type  of  short  range  order  predominates. 
Thus  this  microscopic  theory  for  the  nucleation  resistance  requires  that  atoms  in  the 
liquid  be  positioned  randomly  within  the  repulsive  constraints.  These  atoms  would  then 
have  to  be  localized  sharply  to  form  a  crystal-melt  interface  with  no  density  deficit. 

Effect  of  alloying  on  crystal  nucleation .  There  have  been  several  studies  of  the 
dependence  of  the  temperature  Tc(x),  at  measurable  nucleation  onset,  on  melt 
composition,  x.  They  generally  indicate  [25]  that  in  the  labile  regime  T0(x)  roughly 
parallels  the  liquidus  Tj (x) .  Thompson  and  Spaepen  [26]  showed  that  this  behavior  can 

be  explained  with  simple  solution  models  by  supposing  that:  (i)  The  scaled  crystal- 
melt  tension,  a,  as  in  Spaepen ' s  model,  is  mainly  topological  in  origin  and  only  weakly 
dependent  on  x.  (ii)  The  crystal  nucleus  is  always  in  contact  with  melt  of  average 
composition;  this  supposition  is  the  analog  of  the  "mean  field  approximation". 

3 .  Laser  Quenching 

a.  Quenching  mechanism. 

Irradiating  a  metallic  surface  with  a  short  laser  pulse  results  in  the  melting  and 
subsequent  quenching  of  a  thin  overlay.  It  is  useful  to  analyze  this  process  as  a 
succession  of  three  stages:  energy  deposition,  melt-in,  and  regrowth.  They  are 
illustrated  in  Figure  1. 
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Fig.  1.  Schematic  illustration  of  the  mechanism  of  pulsed  laser  quenching. 

During  the  first  stage,  the  laser  energy  is  deposited  in  a  layer  of  thickness  a-1, 
where  a  is  the  absorption  coefficient.  For  optical  light  on  metals,  or  for  UV 
light  on  semiconductors,  this  layer  is  a  few  hunderd  A  thick.  Since  the  transfer 
of  the  energy  from  the  electrons,  which  interact  with  the  light,  to  the  lattice 
occurs  within  less  than  1  ps  [27],  the  energy  deposition  process,  which  lasts  for 
the  duration  of  the  pulse,  tp,  can  be  considered  a  purely  thermal  one.  The  thermal 

diffusion  length  during  that  time  is  =  {2Dtfttp) ,  which,  for  metals  is 
between  250  and  700  A,  depending  on  the  thermal  diffusivity  (Dth=  10~4  to  10"^ 
m2/s) .  For  fluences,  <}>,  in  the  indicated  range,  a  layer  of  thickness  1T  is  melted 
and  heated  to  several  thousand  degrees. 


Note  that  for  a  pulse  duration  of  less  than  1  ps,  is  less  than  a-^ .  The 
thickness  of  the  initial  molten  layer,  and  hence  also  the  thermal  gradients,  are 
then  independent  of  the  pulse  duration.  Femtosecond  pulses,  therefore,  although 
very  useful  for  time-resolved  spectroscopy  of  the  various  phenomena  associated  with 
the  irradiation  process,  do  not  produce  a  higher  quench  rate.  A  general  discussion 
of  the  regimes  associated  with  the  relative  values  of  lij  and  a-^  can  be  found  in 
Bloembergen ' s  paper  [28]. 

Following  the  termination  of  the  irradiation  pulse,  the  heat  from  the  very  hot 
liquid  continues  to  diffuse  into  the  substrate,  which  results  in  further  melting, 
until  the  temperature  of  the  crystal-liquid  interface  drops  to  the  equilibrium 
melting  temperature.  For  the  conditions  of  Figure  1,  the  maximum  melt  depth,  d,  is 
about  1000A.  The  overheating  of  the  crystal  during  the  melt-in  process  is 

considerable.  The  thermal  gradients  are  very  high  (1011  K/m) ,  as  is  the  melt-in 
velocity  (u-lO^m/s) .  The  duration  of  the  melt-in  stage  of  the  process  is  therefore 
on  the  order  of  d/u  -  100  ps . 


The  third  stage  is  that  of  resolidification,  either  by  regrowth  of  the  underlying 
crystal,  by  growth  of  a  new  crystal,  or  by  glass  formation  if  the  melt  cools  to  its 


configurational  freezing  point  before  the  crystallization  front  passes.  A  number 
of  sophisticated  numerical  analyses  have  been  made  of  the  thermal  quantities  of 
this  process,  but  their  order  of  magnitude  can  easily  be  estimated  from  dimensional 
arguments.  The  temperature  scale  is  set  by  the  melting  temperature,  Tm,  which  is 

on  the  order  of  103K.  The  length  scale  is  set  by  the  melt  depth  d-10-7m.  The 
corresponding  temperature  gradient  is  then  Vt  -  Tm/d  -  lO^K/m,  and  the  cooling 

rate  f  =  Dtjj  VT/d  ~  10*2  K/s.  The  lifetime  of  the  melt  can  be  estimated  as  t  ■ 

Tm/f  =  10-9  s,  which  is  in  agreement  with  the  experimental  values  obtained  from 
transient  reflectivity  measurements  [29-32] .  The  velocity  with  which  the  isotherms 
move  toward  the  surface  during  cooling  is  an  important  quantity  to  be  compared  to 
the  crystal  growth  velocity,  and  can  be  estimated  as  uT  «  T/Vt  -  100  m/s.  These 

parameters,  together  with  those  for  a  "conventional"  melt  quenching  process  such  as 
melt  spinning,  are  summarized  in  Table  II,  using  the  thermal  data  for  iron. 

Table  II.  Thermal  parameters  in  melt  quenching. 


Laser  Quenching  Melt  Spinning 


Melt  temperature 

Tm(K) 

103 

103 

Melt  thickness 

d  (m) 

O 

1 

SxlO-5 

Temperature  gradient 

VT (K/m) 

io*o 

2xl07 

Cooling  rate 

f (K/s) 

10*2 

4xl06 

Melt  lifetime 

t  (s) 

10~9 

not  applicable 

Isotherm  velocity 

Heat-flow  limited 

uij  (m/s) 

100 

0.2 

crystal  growth  velocity 

Uh (m/s) 

230 

0.5 

b.  Thermodynamics  and  kinetics  Qf_ phase  formation- 

The  short  timescale  of  the  picosecond  melt  quenching  process  puts  restrictions  on 
the  transformations  that  can  occur.  For  the  entire  liquid  layer,  of  thickness 
d  «  1000  A,  to  crystallize  during  the  cooling  time,  T,  of  1  ns,  the  velocity  must 
be  at  least  u  =  d/T  -  100  m/s.  The  time,  tj_,  to  crystallize  a  monolayer  at  this 
speed  is  about  3  ps,  and  the  distance  an  atom  can  diffuse  in  that  time,  (Djit^)  *//2, 
is  less  than  an  interatomic  distance  (Dj  is  the  diffusivity  in  the  liquid  phase) . 

This  means  that  under  these  short  pulsed  laser  quenching  conditions  no  long-range 
diffusion  can  occur,  and  that  the  only  possible  transformations  are  partitionless 
ones,  either  growth  of  a  crystal  of  the  same  composition  of  the  liquid,  or  glass 
formation . 

The  thermodynamic  restrictions  imposed  by  the  partitionless  conditions  are 
illustrated  by  the  phase  diagram  of  Figure  2,  which  contains  two  primary  solid 
solutions  (0,8)  and  two  intermetallic  compounds  (|J,y)  .  The  free  energy  diagram  on 
Figure  2  shows  how  the  Tc-lines  for  the  transformations  of  these  four  phases  from 
the  liquid  are  constructed.  At  a  sufficiently  low  temperature,  Tg,  the  atomic 

transport  rate  in  the  liquid  state  becomes  negligibly  small,  so  that  no  more 
transformations  can  occur  on  a  reasonable  timescale.  The  liquid  is  then 
configurationally  frozen  and  becomes  a  glass.  For  a  composition  in  the  range  x'-x" 
in  Figure  2,  no  crystallization  is  possible  under  picosecond  laser  quenching 
conditions,  and  therefore  only  glass  formation  is  thermodynamically  possible.  This 
criterion  for  glass  formation  is  only  a  sufficient  one,  and  not  a  necessary  one,  as 


has  been  claimed  (33] .  In  fact,  more  often  than  not,  picosecond  quenching 
experiments  yield  glasses  over  a  wide  range  of  compositions  below  the  T0-lines  for 

both  intermetallic  compounds  and  primary  solutions.  To  explain  this,  it  is 
necessary  to  consider  the  kinetics  of  partitionless  crystallization  in  some  detail. 
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Fig.  2.  Schematic  phase  diagram  (top),  and  corresponding  free  energy  diagram  at 
temperature  Tj  (bottom),  illustrating  the  construction  of  the  TQ-lines  for  the 
primary  solid  solutions  and  the  intermetallic  compounds.  In  the  composition  range 
x'-x"  (dashed) ,  only  glass  formation  is  possible  under  conditions  of  partitionless 
solidification .  The  numbers  correspond  to  the  same  compositions  in  the  top  and 

bottom  diagrams . 

The  basic  theory  for  crystal  growth  far  from  equilibrium  has  been  formulated  by 
Turnbull,  and  can  be  found  in  a  number  of  publications  [34,35],  including  the  first 
paper  in  the  present  volume.  The  crystal  growth  velocity  can  be  written  as: 

u  -  fkX  [1  -  exp ( - AGc — )]  (15) 


where  k  is  the  atom  jump  frequency  across  the  crystal-melt  interface,  f  is  the 
fraction  of  interface  sites  that  can  incorporate  a  new  atom,  X  is  the  interatomic 
distance,  AGC  is  the  difference  in  molar  free  energy  between  crystal  and  melt,  and 
T i  is  the  interface  temperature. 


If  the  undercooling  of  the  interface,  AT^ 
energy  can  be  approximated  by: 


Tm-Ti,  is  small,  the  driving  free 


AS„  AT; 


where  ASC  is  the  molar  entropy  of  crystallization.  The  exponential  in  equation 
(15)  can  be  linearized  to: 


u  -  fkX  -ASc  AXi 

R  T , 
1 1 


(17) 


where,  for  metals,  f  ~  1  and  ASC  ~  R. 

The  crystal  growth  rate  also  depends  on  the  rate  at  which  the  heat  of 
crystallization,  AhCfV  per  unit  volume,  can  be  removed.  Since  this  must  occur  by 
heat  flow  down  the  temperature  gradient  at  the  interface,  V^T,  the  heat  flux,  J, 
must  obey  the  relations: 

J  =  u  AhC/V  =  K  VjT  (18) 

where  K  is  the  thermal  conductivity.  Conversely,  equation  (18)  can  be  turned  around 
to  define  a  heat  flow-limited  velocity: 

uh  =  -KXiX  (19) 

^c,  v 

Table  I  gives  values  for  u^  for  the  laser  quenching  conditions  as  well  as  for  melt 
spinning. 

Combining  equations  (17)  and  (19)  then  gives  for  metals: 

uh  -  Xm  ~ -Xi  (20) 

kX  Tj. 

In  conventional  solidification,  the  interface  kinetics  are  much  faster  than  the 
rate  of  heat  removal,  so  that  u^  «  kX,  and  AT^  is  small.  This  type  of  growth  is 
called  heat  flow- limited.  In  laser  quenching,  depending  on  the  nature  of  the 
interface  rearrangements  (see  below) ,  the  rate  of  heat  removal  can  be  much  faster 
than  the  interface  kinetics,  so  that  ujj  »  kX,  and  the  undercooling  AT^  is  large, 

so  that  metastable  phases  can  be  formed.  This  corresponds  to  interface-limited 
growth . 

The  maximum  growth  velocity  of  a  crystal,  u^^  -  kX,  is  determined  by  the  nature  of 

the  atomic  rearrangements  that  control  the  growth  process.  In  pure  metals,  dilute 
alloys  or  compounds  with  simple  crystal  structures,  k  can  be  taken  as  the  frequency 
of  thermal  vibration  (~  Debye  frequency),  so  that  u^*  -  u3,  the  speed  of  sound  in 

the  liquid,  which  is  on  the  order  of  several  thousand  m/s.  This  is  called  the 
collision-controlled  regime,  because  the  growth  occurs  by  simple  impingement  of  the 
liquid  atoms  onto  the  crystal. 

In  more  concentrated  alloys  or  in  compounds  with  a  more  complicated  crystal 
structure,  where  a  change  of  the  nearest  neighbor  configuration,  i.e.  diffusion,  is 
necessary,  the  jump  frequency  can  be  taken  as  Dj/X^,  and  =  Dj/X,  which  is  on 

the  order  of  10  m/s.  This  is  the  diffusion-controlled  regime.  Table  I  shows  that 
for  picosecond  laser  quenching  umax>D  <<  uT,  u^,  so  that  this  type  of  growth  can 
easily  be  suppressed.  Equation  (20)  shows  that  under  these  conditions  T^  falls  far 
below  Tm,  and  glass  formation  becomes  possible. 


Although  in  partitionless  growth  of  crystals  there  is  no  long-range  diffusional 
transport,  the  growth  can  still  be  diffusion-controlled  if  a  drastic  change  of  the 
nearest  neighbor  environment  is  necessary  to  transform  the  liquid  structure  into 
the  crystalline  one.  This  is  the  case,  for  example  in  metal-metalloid  alloys  such 
as  Fe-B,  where  the  chemical  interaction  leads  to  a  very  pronounced  short  range 
order  in  the  liquid  <9  Fe  atoms  around  each  B  in  a  capped  trigonal  prism 
arrangement  [36]),  which  is  very  different  from  that  in  the  crystal  (B  interstitial 
in  bcc  Fe  [37]).  As  a  result,  it  is  possible  to  produce  Fe-B  glasses  containing  as 
little  as  5  at . %  B  [38],  which  is  far  below  the  T0-line.  Similarly,  if  the  unit 

cell  of  an  intermetallic  compound  is  large,  or  if  long-range  chemical  order  must  be 
established,  crystal  growth  is  diffusion-limited,  and  the  growth  of  these  compounds 
is  easily  suppressed  in  picosecond  pulsed  laser  quenching.  Numerous  experiments  on 
a  variety  of  binary  metallic  systems  (Ni-Nb,  Ni-Mo,  Co-Mo,  Nb-Co  [30,39-43])  have 
demonstrated  that  glasses  can  be  formed  at  the  intermetallic  compositions.  On  the 
other  hand,  in  pure  metals  with  simple  crystal  structures,  and,  by  extension, 
disordered  alloys  with  the  same  simple  structures  (e.g.  disordered  fee  in  Cu-Co  or 
Ag-Cu) ,  growth  is  collision-limited  and  cannot  be  suppressed  even  by  picosecond 
pulsed  laser  quenching.  This  is  confirmed  by  direct  measurements  of  the  growth 
velocity  of  pure  metals  and  disordered  alloys  (Au,  Cu,  Cu-Au)  by  transient 
reflectivity  [30,31]  which  yield  large  values  (>100  m/s)  that  can  only  be  accounted 
for  by  a  collision-controlled  mechanism. 
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Fig.  3.  Phase  diagram  of  the  Nb-Si  system  (top),  and  schematic  illustration  of  the 
phases  obtained  by  different  methods  of  liquid  quenching.  The  melt  spinning 
results  are  from  refs.  [48-51];  the  laser  quenching  results  are  from  refs.  [45-4 7], 
The  pulse  length  and  the  substrate  for  the  laser  experiments  are  indicated. 


An  undercooled  liquid  can  often  crystallize  into  any  of  a  number  of 
thermodynamically  allowed  stable  or  metasable  phases.  The  selection  of  the  phase 


that  actually  forms  depends  on  the  details  of  the  transformation  kinetics.  Some 
general  empirical  rules  to  predict  this  have  been  formulated.  The  best  known  one 

of  these  is  Ostwald's  rule,  which  states  that  the  phase  nearest  to  the  liquid  in 

free  energy  is  kinetically  favored.  There  are,  however,  numerous  counterexamples 
to  this  rule.  Turnbull  has  proposed  that  phases  closest  in  entropy  to  the  parent 
phase  are  kinetically  favored  (44]  .  He  pointed  out  that  that  transformation 
processes  requiring  the  smallest  changes  in  the  positional  and  maitonal 
correlations  of  the  atoms  are  likely  to  be  the  easiest  ones.  For  tranformations 
from  the  liquid  this  implies  that  a  small  unit  cell  and  a  high  degree  of  isotropy 
are  favorable  kinetic  factors  for  crystallization. 

A  recent  pulsed  laser  quenching  study  of  the  Nb-Si  system  illustrates  these  factors 
very  clearly  [45-47].  The  results  are  summarized  in  Figure  3. 

At  the  highest  quench  rates  (>10^  k/s) ,  obtained  in  the  thinnest  sections  of  the 

spots  irradiated  with  a  30  ps  pulse  on  a  copper  substrate  the  amorphous  phase  was 
formed  over  a  wide  composition  range.  Such  broadening  of  the  range  of  amorphous 
phase  formation  is  consistent  with  our  earlier  observations  on  other  metal- 
metalloid  and  metal-metal  systems. 

At  slightly  lower  quench  rates  (10^2  K/s),  obtained  in  the  remaining  area  of  those 
spots,  a  new  phase  was  found:  a  disordered  f.c.c.  solid  solution  [46].  Although 
the  formation  of  disordered  f.c.c.  solid  solutions  in  metal-metal  alloys  cannot  be 
suppressed  by  picosecond  laser  quenching,  the  growth  rate  of  the  f.c.c.  phase  in 
this  system,  although  still  the  highest  one  observed,  is  probably  slightly  lower 
due  to  the  short  range  order  around  the  Si. 

By  lowering  the  quench  rate  to  lO^-lO1^  k/s  by  supporting  the  films  on  an  aluminum 
substrate,  a  supersaturated  b.c.c  solid  solution  was  found  at  low  Si  contents,  and 
an  ordered  f.c.c.  phase  (LI2,  Cu3Au-type)  around  25  at.%  Si. 

A  cooling  rate  of  109-1010  K/s  was  obtained,  either  by  supporting  the  ps-irradiated 
films  on  a  salt  substrate,  or  by  irradiating  the  films  on  a  copper  substrate  with  a 
30  ns  pulse.  Under  these  conditions  two  entirely  new  phases  were  found  around  25 
at.%  Si:  a  hexagonal  phase  with  a  250  unit  cell  and  an  as  yet  unidentified  phase 
[47], 

Conventional  melt  spinning,  during  which  the  cooling  rate  is  10®  K/s,  produces  the 
TijP-type  structure  [48], 

Important  to  note  here  is  that  the  degree  of  isotropy  and  the  unit  cell  size  of  the 
phases  formed  at  25  at%  Si  decrease  with  decreasing  quench  rate:  amorphous 
(isotropic,  no  unit  cell),  disordered  fee  (1  atom  unit  cell,  cubic  point  group), 
ordered  f.c.c  (4  atom  unit  cell,  cubic  point  group),  hexagonal  (-  16  atom  unit 
cell,  uniaxial  point  group),  TijP-type  (32  atom  unit  cell,  uniaxial  point  group). 

Qualitatively  similar  observations  have  been  made  in  laser  quenching  of  Ga  [52]  and 
Mn  [53]  . 
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Abstract 

Laser  quenching  was  used  to  investigate  the  forma¬ 
tion  of  metastable  phases  in  the  Nb-Si  allay  system  in 
the  range  8-27  at.  %  Si.  By  varying  the  duration  of  the 
laser  pulse  between  3(J  ps  and  30  ns.  and  by  changing 
the  conductivity  of  the  substrate,  melt  quench  rates  be¬ 
tween  10 8  K  s  1  and  10 K  s  1  were  obtained.  At  the 
highest  cooling  rates,  amorphous  alloys  were  formed 
over  the  entire  composition  range.  At  progressively 
lower  quench  rates  the  following  metastable  phases  are 
/armed:  face -centered  cubic  solutions,  supersaturated 
body -centered  cubic  solutions,  the  LI,  ordered  phase, 
and  a  new  hexagonal  compound.  This  progression  illus¬ 
trates  how  phase  formation  from  the  melt  is  kineticalh 
favored  by  higher  isotropy  and  a  smaller  unit  cell  size. 
All  the  phases  exhibited  superconductivity,  with  a  T 
between  2  K  and  6  K. 

I.  Introduction 

Pulsed  laser  quenching  is  the  fastest  melt  quench¬ 
ing  technique  available:  picosecond  pulses  can  give 
quench  rates  of  the  order  of  10i:  K  s  The  technique 
has  been  used  to  investigate  the  compositional  limits 
of  glass  formation  and  to  expose  new  metastable 
phases  [1-4].  Until  very  recently,  the  stoichiometric 
A- 15  phase  of  Nb-Si  was  the  prime  theoretical  candi¬ 
date  for  the  highest  critical  temperature.  Since  this 
phase  is  metastable.  it  must  be  prepared  by  non-equi¬ 
librium  methods.  Conventional  melt  spinning  can 
produce  the  A- 15  structure,  but  only  if  it  is  niobium 
rich  At  25  at  °o  Si.  it  produces  the  tetragonal  Ti,P 
structure  [5-8], 

Although  the  laser  quenching  experiments  de¬ 
scribed  here  did  not  produce  the  A- 1 5  phase,  they  led 
to  the  discovery  of  some  new  metastable  phases  and 
provided  insight  into  the  fundamentals  of  metastable 
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phase  selection.  This  paper  explains  the  basics  of  laser 
quenching  and  phase  selection,  and  reports  the  results 
of  its  application  to  Nb-Si. 


2.  Laser  quenching 

Irradiating  a  metallic  surface  with  a  short  laser 
pulse  results  in  the  melting  and  subsequent  quenching 
of  a  thin  overlay.  This  process  can  be  analyzed  as  a 
succession  of  three  stages:  energy  deposition,  melt-in 
and  regrowth  [1.9].  They  are  illustrated  in  Fig.  1  for 
a  pulse  of  30  ps. 

During  the  first  stage,  the  laser  energy  is  deposited 
as  heat  in  a  layer  of  thickness  a  where  i  is  the 
absorption  coefficient,  which  for  metals  is  a  few  hun¬ 
dred  angstroms.  The  thermal  diffusion  length  during 
that  time  is  / T  =■  (2£>lh/p)'  which  is  between  250  A 
and  700  A  for  metals,  depending  on  the  thermal  diffu- 
sivity.  For  fluences  <p  in  the  indicated  range,  a  layer  of 
thickness  lr  is  melted  and  heated  to  several  thousand 
degrees. 

After  the  pulse,  the  heat  from  the  very  hot  liquid 
continues  to  diffuse  into  the  substrate,  resulting  in 
further  melting,  until  the  temperature  of  the  crystal- 
liquid  interface  drops  to  the  equilibrium  melting  tem¬ 
perature.  For  the  conditions  of  Fig.  1.  the  maximum 
melt  depth  d  is  about  1000  A. 

The  third  stage  is  that  of  resolidification,  either  by 
regrowth  of  the  underlying  crystal,  by  growth  of  a 
new  crystal,  or  by  glass  formation  if  the  melt  cools  to 
its  configurational  freezing  point  before  the  crystal¬ 
lization  front  passes.  The  quench  parameters  can  be 
estimated  by  the  following  dimensional  argument. 
The  characteristic  temperature  is  the  melting  tem¬ 
perature  Tm  which  is  of  the  order  of  10’  K.  The  char¬ 
acteristic  length  is  the  melt  depth  d  =  10  ’  m. 
The  corresponding  temperature  gradient  is  then 
NT  =  Tm  d  i  10'°  K  m  '.  and  the  cooling  rate 
t  =  Dth  NT  d  ~  10'-  Ks  '.  The  lifetime  of  the  melt  is 
t  —  Tm  t  =  10  gs.  which  is  also  the  value  found  in 
transient  reflectivity  measurements  [  10,  1 1], 
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I  Schematic  illustration  of  the  mechanism  of  pulsed  laser  quenching 


3.  Phase  selection 

The  short  timescale  of  the  picosecond  melt  quench¬ 
ing  process  puts  restrictions  on  the  transformations 
that  can  occur.  For  the  entire  liquid  layer,  of' thickness 
</  =  1000  A.  to  crystallize  during  the  cooling  time  t  of 
I  ns.  the  velocity  must  he  at  least  u  =  d  i  ±  100  ms  '. 
The  time  r,  to  <-r;.  stallize  a  monolayer  at  this  speed  is 
about  .1  ps.  and  the  distance  an  atom  can  diffuse  in 
that  time.  {D,t ,  I '  \  is  less  than  an  interatomic  distance 
( D,  is  the  diffusivity  in  the  liquid  phase).  This  means 
that  under  these  short  pulsed  laser  quenching  condi¬ 
tions.  no  long-range  diffusion  can  occur  and  the  only 
possible  transformations  are  partitumless  ones,  either 
growth  of  a  crystal  of  the  same  composition  of  the 
liquid,  or  glass  formation.  For  nanosecond  pulses  or 
for  high  fluences.  which  result  in  a  larger  melt  depth, 
the  regrowth  velocity  is  lower,  and  some  degree 
ol  phase  separation,  usually  on  a  very  fine  scale, 
become^  possible 

In  conventional  solidification,  the  interface  kinet¬ 
ics  are  much  faster  than  the  rate  of  heat  removal,  so 
that  the  undercooling  at  the  interface.  AT,,  is  small. 
This  type  of  growth  is  called  heat  How  -limited  In  laser 
quenching,  depending  on  the  nature  of  the  interlace 
rearrangements,  the  rate  of  the  heat  removal  can  be 
much  faster  than  the  interface  kinetics,  so  that  AT,  is 
large,  and  meiastable  phases  can  be  formed  This  cor¬ 
responds  to  interface-limited  growth  [1.  12.  I3j 

The  number  of  metastable  phases  with  a  free  en¬ 
ergy  lower  than  that  of  the  liquid  can  be  considerable 


at  large  undercooling,  even  under  the  constraint  of 
partitionless  growth.  The  principle  that  governs  the 
selection  among  them  has  been  the  subject  of  continu¬ 
ing  debate.  The  best  known  of  these  is  the  Ostwald 
step  rule,  which  states  that  the  transformation  to  the 
phase  closest  in  free  energy  is  most  likely  to  occur. 
Turnbull  has  recently  suggested  that  it  is  the  phase 
that  is  closest  in  entropy  that  is  kinetically  favored 
[14],  Fie  pointed  out  that  transformation  processes 
that  require  the  smallest  changes  in  the  positional  and 
motional  correlations  are  likely  to  be  the  easiest  ones. 
For  a  crystallization  process,  this  means  that  phases 
with  a  small  unit  cell  and  a  high  degree  of  isotropy  are 
favored  to  nucleate  or  grow 

4.  Experiments 

The  starting  samples  for  laser  irradiation  were  pre¬ 
pared  as  artificial  multilayers  by  sputtering  alternat¬ 
ing  targets  |4.  I  5).  The  substrates  were  either  copper, 
aluminum  or  KCI.  The  geometry  is  shown  schemati¬ 
cally  in  Fig.  2.  The  periodicity  of  the  layers  was  kept 
smaller  than  3  nm.  which  is  the  difTusional  mixing 
length  in  the  liquid  for  the  I  ns  lifetime  of  the  melt. 
The  individual  layers  were  elemental  niobium  (poly- 
crystalline  body-centered  cubic  ( b.c.c. ) )  and  elemental 
silicon  (amorphous),  for  the  low  silicon  concentra¬ 
tions.  the  silicon  target  was  replaced  by  a  composite 
Nb  Si  target  The  average  composition  of  the  samples 
could  be  adjusted  by  changing  the  relative  thicknesses 
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Fig.  2  Schematic  diagram  (on  a  variable  scale)  of  a  multilayer 
film  after  irradiation  with  a  laser  pulse  with  a  gaussian  intensity 
profile 


of  the  layers,  and  covered  the  range  8-27  at.%  Si. 
After  irradiation  with  a  30  ps  Nd:YAG  laser  or  with  a 
30  ns  ruby  laser,  the  films  were  chemically  removed 
Irom  the  substrate.  They  could  then  be  examined  in 
the  transmission  electron  microscope  without  further 
thinning. 

The  superconducting  transition  temperature  7 1 
was  determined  on  samples  that  had  a  line  of  irradia¬ 
tion  spots  between  the  electrodes.  The  transition 
could  be  identified  by  a  drop  in  the  resistance. 
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fig  3  Phase  diagram  of  the  Nh  Si  system  ( top),  and  schematic 
illustration  of  the  phases  obtained  by  difTercnt  methods  of  liq¬ 
uid  quenching  The  melt-spinning  results  are  from  refs  5  8. 


5.  Results 

The  results  of  the  quenching  experiments  are  sum¬ 
marized  in  Fig.  3.  At  the  highest  quench  rates  (above 
10'-  Ks  1 ).  obtained  in  the  thinnest  sections  of  the 
spots  irradiated  with  a  30  ps  pulse  on  a  copper  sub¬ 
strate.  the  amorphous  phase  was  formed  over  the  en¬ 
tire  composition  range  investigated.  This  is  a  much 
wider  range  of  amorphous  phase  formation  than  is 
possible  '  y  other  melt-quenching  methods.  This  ob¬ 
servation  is  consistent  with  our  earlier  observations 
on  other  metal-metalloid  [2]  and  metal  metal  [4] 
systems. 

At  slightly  lower  quench  rates  (  1 0 1  -  K  s  ').  ob¬ 
tained  in  the  remaining  area  of  those  spots,  a  new 
phase  was  found:  a  disordered  fc.c.  solid  solution 
[  16]  Although  the  formation  of  disordered  f.c.c.  solid 
solutions  in  metal  metal  alloys  cannot  be  suppressed 
by  picosecond  laser  quenching  [  I,  3],  the  growth  rate 
of  the  f.c.c.  phase  in  this  system,  although  still  the 
highest  observed,  is  probably  slightly  reduced  because 
of  the  short-range  order  around  the  silicon  [2j. 

By  lowering  the  quench  rate  to  ( 10"  1  O'  * )  K  s  ' 
by  supporting  the  films  on  an  aluminum  substrate,  a 
supersaturated  b.c.c.  solid  solution  was  found  at  low 
silicon  contents,  and  an  ordered  f.c.c.  phase  (LI,. 
Cu,Au  type)  at  around  25  at.%  Si 


A  cooling  rate  of  ( I  O'*  10"’)  K  s  1  was  obtained, 
either  by  supporting  the  picosecond-irradiated  films 
on  a  salt  substrate,  or  by  irradiating  the  films  on  a 
copper  substrate  with  a  30  ns  pulse.  Under  these  con¬ 
ditions.  two  entirely  new  phases  were  found  around 
25  at .%  Si:  a  hexagonal  phase  with  a  250  A'  unit  cell 
and  an  as-vet  unidentified  phase  [  17]. 

7  _ _ , _ _ _ _ _ _ _ _ _ _ _ _ _ , 


Fig  4  Superconducting  transition  temperature  as  a  function  of 
silicon  content  for  the  as-prepared  modulated  Nb  Si  films,  and 
for  the  amorphous,  bcc  and  fee  phases 
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The  results  of  the  superconductivity  measurements 
are  summarized  in  Fig.  4.  A  general  decrease  of  Tc 
with  silicon  content  is  observed.  Tc  values  for  the  new 
hexagonal  phase  and  the  as-yet  unidentified  phases 
were  2.7  K  and  2.3  K  respectively. 

6.  Discussion 

As  an  illustration  of  the  phase  selection  rules  dis¬ 
cussed  above,  it  is  interesting  to  note  the  degree  of 
isotropy  and  the  unit  cel!  size  of  the  phases  formed  at 
25  at.%  Si  in  order  of  decreasing  quench  rate:  amor¬ 
phous  (isotropic,  no  unit  cell),  disordered  f.c.c.  (one- 
atom  unit  cell,  cubic  point  group),  ordered  f.c.c. 
(four-atom  unit  cell,  cubic  point  group),  hexagonal 
(unit  cell  of  about  16  atoms,  uniaxial  point  group). 
Ti,P  type  (32-atom  unit  cell,  uniaxial  point  group). 
Increasing  the  quench  rate  increases  the  undercooling 
and  exposes  phases  that  are  increasingly  more 
metastable.  It  is  clear  from  the  experiments  that  the 
more  isotropic  phases  with  a  smaller  unit  cell  are 
kinetically  favored. 

The  decrease  of  7"v  with  silicon  content  in  the  mod¬ 
ulated  films  was  also  observed  by  DenhofT  and  Gygax 
[18].  A  similar  decrease  in  the  amorphous  phase  may 
be  the  result  of  the  increase  in  the  number  of  electrons 
per  atom,  as  proposed  by  Collver  and  Hammond  for 
metallic  glasses  ( 19).  The  same  explanation  may  hold 
for  the  f.c.c.  and  b.c.c.  solid  solutions.  The  value  of 
for  the  f.c.c.  solution  found  here  is  almost  the  same  as 
that  found  for  the  same  phase  produced  by  shock 
compression  [20]  It  is  interesting  to  note  that  Tt  for 
the  metastable  f.c.c.  solutions  is  higher  than  that  for 
the  b.c.c.  solutions  by  roughly  the  same  amount  over 
their  entire  composition  ranges.  One  might  speculate 
that  T  for  metastable  pure  f.c.c.  niobium  should  be 
higher  than  the  9.25  K  value  for  the  stable  b.c.c. 
phase.  The  TK  values  for  the  hexagonal  and  uniden¬ 
tified  phases  are  relatively  high,  since  the  non-cubic 
Nb  Si  compounds  usually  have  7",  values  below  I  K. 

Acknowledgments 

We  thank  J.  L.  Bell  for  assistance  with  the  sample 
preparation,  and  A  M.  Malvezzi  and  C.  A.  MacDon¬ 


ald  for  help  with  the  laser  irradiation.  This  work  has 
been  supported  by  the  Office  of  Naval  Research  under 
contract  NOOOI4-K-85-0684.  W.K.W.  has  been  a 
participant  in  the  Chinese-American  Cooperative 
Basic  Research  Program  in  Atomic,  Molecular  and 
Condensed  Matter  Physics  of  the  American  Physical 
Society. 


References 

1  F.  Spaepen.  in  E  W  Collings  and  C.  C.  Koch  (eds). 
Undercooled  Alloy  Phases ,  TMS-AIME.  New  York,  1987, 
p.  187. 

2  C.  J  Lm  and  F.  Spaepen.  Appl  Phys.  Lett.,  41  (1982) 

721. 

3  C.  J  Lin.  F  Spaepen  and  D  Turnbull.  J  Non-Crvsl. 
Solids.  61  62  ( 1984)  767 

4  C  J  Lin  and  F.  Spaepen.  Acta  Meiall .  34(1986)  1367 

5  A  lnoue  and  T.  Masumolo.  Sci.  Rep  Res.  Inst  Tohoku 
University.  A.  28(1980)  165. 

6  R  M  Watersirat,  F.  Haenssler  and  J  Muller.  J  Appl. 
Phys  .  SO  ( 1979) 4763 

7  K  Togano.  H.  Kamakura  and  K  Tachikawa.  Phys  Lett. 
A.  76(1980)  83. 

8  W  Rossteutscher  and  K  Schubert.  7  Metallkd..  56(1965) 
813 

9  N  Bloembergen.  in  S  D  Ferns.  H  J  Leamv  and  J  M 
Poate  (eds  ).  Laser-Solid  Interactions  and  Laser  Processing. 
A1P.  New  York.  1979.  p  i 

10  J  M  Liu.  R  Yen.  H.  Kurz.  and  N  Bloembergen.  Appl 
Phis  Lett..  39(1981)  755 

1 1  C  A  MacDonald.  A  M  Malvezzi  and  F  Spaepen.  Mater. 
Res  Six  Symp  Proc..  SI  11986)  277. 

12  D  Turnbull,  m  M  Beriollolti  (ed  ).  Physical  Processes  in 
Laser-Materials  Interactions,  NATO-AS1  Proc..  Plenum. 
New  York.  1983,  p  117. 

13  F  Spaepen  and  D.  Turnbull,  in  J  M  Poate  and  J.  W. 
Mayer  (eds).  Laser  Processing  of  Semiconductors. 
Academic  Press.  New  York.  1982.  p  15 

14  D  Turnbull.  Meiall.  Trans  A.  /2(I98I)  695 

15  F  Spaepen,  A  L  Greer.  K  F  Kelton  and  J  L  Bell.  Ret 
Sci  Instrum  .  561  1985)  1343. 

16  W  K  Wang  and  F  Spaepen.  J  Appl  Phys..  Sfi  ( 1985) 
4477. 

1 7  W  K  Wang  and  F.  Spaepen.  Appl  Phys.  Lett..  49  ( 1986) 
295 

18  M  DenhofT  and  G  Gygax.  in  U  Eckern.  A  Schmidt.  W 
Weber  and  H  Wuhl  (eds  ).  Proc  Conf  on  Lon  Tempera¬ 
tures.  LTI7.  Elsevier.  Amsterdam.  1984.  p  1317 

19  M  M  Collver  and  R  H  Hammond.  Phys.  Rex  Lett  .  30 
( 1973)  92 

20  W  K  Wang.  Y  Svono.  G  Goto.  H  Iwasaki.  A  lnoue 
and  T  Masumolo.  Scr.  Meiall .  /5(|98I)  1313 


